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EXECUTIVE  SUMMARY 


As  part  of  the  overall  NCEL  marine  cable  dynamics  exploratory  development  program  under  the 
sponsorship  of  NAVFAC,  a  series  of  laboratory  and  field  experiments  have  been  conducted  to  investi¬ 
gate  the  effects  of  attached  masses  and  sensor  housings  (discrete  or  lumped  masses)  on  the  overall 
cable  system  response.  Towing  channel  experiments  were  conducted  with  a  "strumming  rig"  developed 
for  the  NAVFAC/NCEL  cable  dynamics  program  and  the  test  findings  recently  were  reported  (3)  ’  A 
field  lest  program  was  conducted  during  the  summer  of  1981  to  investigate  further  the  strumming 
vibrations  of  marine  cables  in  a  controlled  environment.  The  experiments  were  funded  by  NCEL,  the 
USGS  and  industry  sponsors,  planned  by  NRL  and  MIT,  and  conducted  at  the  field  site  by  MIT  The 
primary  objective  of  the  test  program  was  to  acquire  data  to  validate  and,  if  necessary,  to  provide  a  basis 
for  modifying  the  NCEL-sponsored  computer  code  NATFREQ  (4)  This  code  was  developed  at  the 
California  Institute  of  Technology  to  provide  a  means  for  calculating  the  natural  frequencies  and  mode 
shapes  of  taut  marine  cables  with  large  numbers  of  attached  masses. 

The  vortcx-cxcitcd  oscillations  of  marine  cables,  commonly  termed  strumming,  result  in  increased 
fatigue,  larger  hydrodynamic  forces  and  amplified  acoustic  flow  noise,  and  sometimes  lead  to  structural 
damage  and  eventually  to  failure.  Flow-excited  oscillations  very  often  are  a  critical  factor  in  the  design 
of  underwater  cable  arrays,  mooring  systems,  drilling  risers,  and  offshore  platforms.  Many  components 
of  these  complex  structures  typically  have  a  cylindrical  cross-section  which  is  conducive  to  vortex  shed¬ 
ding  when  they  are  placed  in  a  flow.  An  understanding  of  the  nature  of  the  fluid-structure  interaction 
which  results  in  vortex-excited  oscillations  is  an  important  factor  in  the  reliable  design  and  long-term 
operation  of  offshore  structures  and  cable  systems. 

As  the  state  of  the  ocean  engineering  art  steadily  progresses,  more  and  more  stringent  demands 
are  being  placed  upon  the  performance  of  cable  structures  and  moorings.  In  particular,  displacement 
tolerances  and  constraints  in  response  to  currents  are  becoming  more  stringent;  fatigue  is  becoming  an 

important  design  consideration;  and  the  sensitivity  of  acoustic  sensors  has  increased  to  the  point  that 
'Numbers  in  parentheses  denote  references  listed  el  the  end  of  this  report 


they  cannot  differentiate  between  legitimate  acoustic  targets  and  slight  variations  in  their  vertical  posi¬ 
tion.  All  of  these  are  problems  that  are  aggravated  by  cable  strumming.  In  order  to  design  a  structure 
or  cable  system  to  meet  the  constraints  imposed  by  operational  and  environmental  requirements,  the 
engineer  must  be  able  to  assess  the  effect  of  strumming  on  the  structure  in  question. 

One  purpose  of  this  report  is  to  describe  the  field  test  program  and  to  present  the  salient  results 
from  it.  Time  histories  of  the  measured  hydrodynamic  drag  coefficients,  current  speeds,  and  cable 
strumming  accelerations  and  displacement  amplitudes  are  presented  and  discussed  to  the  extent  possi¬ 
ble.  The  natural  frequencies  and  mode  shapes  for  the  cable  with  attached  masses  were  measured  in  air 
and  in  water  Data  from  the  twenty  test  runs  which  were  conducted  are  presented  in  this  report. 

Validation  of  the  NATFREQ  code  is  an  objective  of  this  report.  Calculations  of  cable  natural  fre¬ 
quencies  and  mode  shapes  using  the  NATFREQ  code  were  made  at  NRL  for  all  of  the  field  test  runs. 
A  comparison  is  made  in  this  report  between  the  calculated  natural  frequencies  and  mode  shapes  and 
those  experimental  test  runs  which  it  was  possible  to  analyze  in  sufficient  detail.  Also,  recommenda¬ 
tions  are  made  of  modifications  that  will  enhance  the  utility  and  the  ease  of  access  to  the  code  for  pros¬ 
pective  users. 

The  test  runs  which  were  conducted  during  the  experimental  phase  of  the  test  program  consisted 
of  ten  pairs  of  equivalent  tests  in  air  and  in  water.  The  measured  in-.;ir  natural  frequencies  are  in  good 
agreement  with  the  NATFREQ  predictions  for  the  second  and  higher  (up  to  n  —  5)  cable  modes  The 
first  mode  frequency  apparently  was  influenced  by  the  sag  of  the  cable.  The  measured  mode  shapes  of 
the  cable  vibrations  in  air  are  in  agreement  with  the  computed  mode  shapes,  but  only  limited  mode 
shape  comparisons  are  possible  due  to  the  existing  capabilities  of  the  code. 

Good  agreement  also  was  obtained  between  the  NATFREQ-prcdicted  natural  frequencies  and  the 

frequencies  measured  in  water  It  is  clear  from  the  comparison  given  in  this  report  that  it  is  possible  to 

predict  the  natural  frequencies  of  a  taut  cable  with  attached  masses  to  better  than  10  percent  However, 

at  high  cable  mode  numbers  the  difference  in  frequency  between  modes  often  is  considerably  less  than 

r. 
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this  value  Then  it  is  difficult  to  positively  identify  which  measured  cable  frequency  is  associated  with 
which  predicted  mode. 

The  results  available  to-date  tend  to  validate  the  NATFREQ  computer  code  as  a  reliable  engineer¬ 
ing  model  for  predicting  the  natural  frequencies  and  mode  shapes  of  taut  marine  cables  with  arrays  of 
discrete  masses  attached  to  them. 

Measurements  of  the  hydrodynamic  drag  forces  during  the  Castine  Bay  field  experiments  con¬ 
sistently  produced  drag  coefficients  in  the  range  Q,  -  2  to  3.2  for  the  bare  cable  and  the  cable  with 
attached  masses.  This  is  a  substantial  amplification  of  the  drag  force  from  the  expected  level  for  a  sta¬ 
tionary  cable  (Cp  -  I  to  I  SI 

The  drag  coefficient  on  the  strumming  cable  was  predicted  by  Vandiver  (23)  with  a  strumming 
drag  model  that  was  developed  at  NRL  as  part  of  the  overall  NCEL  cable  dynamics  research  program 
The  predicted  drag  on  the  cable  was  comparable  to  the  measured  drag.  For  a  complementary  experi¬ 
ment  at  the  test  site  with  a  flexible  circular  steel  pipe,  the  predicted  drag  coefficients  were  virtually 
indistinguishable  from  the  measured  drag  coefficients  (23). 

Several  recommendations  have  been  developed  as  a  result  of  the  comparison  between  the 
NATFREQ  code  predictions  and  the  Castine  Bay  field  tests.  The  major  recommendations  are: 

•  Cable  sag  or  slack  effects  often  play  an  important  role  in  the  dynamics  of  marine  cables 
The  NATFREQ  code  is  limited  to  taut  cable  dynamics.  Consideration  should  be  given  to 
the  development  of  a  comparable  code  for  the  calculation  of  the  natural  frequencies  and 
mode  shapes  of  a  slack  cable  with  attached  masses.  At  least  one  such  code,  called  SLAK. 
presently  exists  in  rudimentary  form  (l). 

•  Many  of  the  data  records  from  the  1981  Castine  field  experiment  contain  lengthy  time 
segments  where  the  cable  strumming  is  nonresonant,  i.e.,  the  vibrations  and  the  vortex 
shedding  are  not  locked-on  at  a  single  resonant  frequency.  Consideration  should  be  given 
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to  developing  a  simple  but  still  effective  method  for  taking  nonresonant  strumming  effects 
into  account  in  determining  the  cable  response  and  the  strumming- induced  hydrodynamic 
drag  force  and  coefficient  CD. 


I.  INTRODUCTION 


1.1  Objectives 

At  the  beginning  of  FY  1975,  technical  coordination  of  the  NAVFAC  marine  cable  dynamics 
exploratory  development  program  was  undertaken  by  the  Naval  Civil  Engineering  Laboratory  (NCEL) 
The  overall  objective  of  this  program,  as  stated  in  a  research  plan  put  forward  by  NCEL,  is: 

"...  10  provide  for  ihe  development  of  effective  methods  for  the  analysis  of  the  dynamic  response 
of  3-dimensional,  moored  cable  structures  which  undergo  dynamic  motions  generated  by  various 
natural  or  man-produced  causes.  Failure  to  predict  this  dynamic  behavior  by  suitable  analytical  tech¬ 
niques  will  affect  the  confidence  in  the  adequacy  of  the  system  design  as  well  as  t he  esnmated  reliability 
of  the  system's  performance  " 

"Dynamic  response  of  moored  cable  structures  includes  two  specific  problem  areas  which  are 
addressed  separately  in  this  plan.  (One  aspect  of  this  problem]  is  the  small-displacement,  "high  fre¬ 
quency"  response  generated  by  shedding  vortices  as  water  flows  past  the  cable  — this  response  is  com¬ 
monly  referred  to  as  cable  strumming.  The  objective  of  the  plan  for  this  specific  area  is  twofold:  (1) 
development  of  a  capability  to  predict  the  strumming  response  (i.e..  deflection,  frequency,  generated 
acoustic  energy,  and  drag  force)  of  cable  networks  which  have  horizontally  or  vertically  oriented  cable 
segments,  m  taut  or  catenary  configurations  subjected  to  a  current  which  may  vary  along  the  cable 
length,  and  (2)  development  of  techniques  which  can  be  used  to  suppress  cable  strumming  " 

The  resonant  strumming  response  of  bare  cables  is  discussed  in  detail  in  a  recent  NCEL/NRL 
report  (I  >  The  suppression  of  strumming  vibrations  is  dealt  with  in  a  separate  NCEL -sponsored  report 
(2) 


As  part  of  the  overall  NCEL  cable  dynamics  program,  both  laboratory  and  field  experiments  have 
been  conducted  to  investigate  the  effects  of  attached  masses  and  sensor  housings  (discrete  or  lumped 
masses)  on  the  overall  cable  system  response.  Towing  channel  experiments  were  conducted  with  a 


'strumming  rig*  developed  for  the  NCEL  cable  dynamics  program  and  the  test  findings  recently  were 
reported  (3).  A  test  program  was  conducted  during  the  summer  of  1981  to  investigate  further  the 
strumming  vibrations  of  marine  cables  in  a  controlled  field  environment.  The  experiments  were  funded 
by  NCEL,  the  USGS  and  industry  sponsors,  planned  by  NRL  and  MIT,  and  conducted  at  the  field  site 
by  MIT.  A  primary  objective  of  the  test  program  was  to  acquire  data  to  validate  and,  if  necessary,  to 
provide  a  basis  for  modifying  the  NCEL-sponsored  computer  code  NATFREQ  (4).  This  code  was 
developed  at  the  California  Institute  of  Technology  to  provide  a  means  for  calculating  the  natural  fre¬ 
quencies  and  mode  shapes  of  taut  marine  cables  with  large  numbers  of  attached  masses. 

1.2  Background 

It  often  is  found  that  bluff,  or  unstrcamlincd.  cylindrical  structures  undergo  damaging  oscillatory 
instabilities  when  wind  or  water  currents  flow  over  them.  A  common  mechanism  for  resonant,  flow- 
excited  oscillations  is  the  organised  and  periodic  shedding  of  vortices  as  the  flow  separates  alternately 
from  opposite  sides  of  a  long  bluff  body.  The  flow  field  exhibits  a  dominant  periodicity  and  the  body  is 
acted  upon  b>  time-varying  pressure  loads.  These  pressures  result  in  steady  and  unstcad)  drag  or  resis¬ 
tance  forces  in  line  with  the  flow  and  unsteady  lift  or  side  forces  normal  to  the  fl'>w  direction  If  the 
structure  is  flexible  and  lightly  damped  as  in  the  case  of  a  marine  cable,  then  .esonant  oscilla’ions  can 
be  excited  normal  or  parallel  to  the  incident  (low  direction  l  or  the  more  common  cross  flow  oscilla¬ 
tions.  the  body  and  the  wake  have  the  same  frequency  of  oscillation  whie  t  is  near  one  ol  the  natural 
frequencies  of  the  system.  The  shedding  meanwhile  is  shifted  away  Iro  n  the  usual,  or  Strmihal.  fre¬ 
quence  at  which  pairs  of  vortices  would  he  shed  if  the  structure  or  cabl  was  restrained  from  oscillating 
This  phenomenon  is  known  as  "lock-in"  or  "wake  capture  " 

The  oscillations  of  marine  cables  caused  by  vortex  shedding,  commonly  termed  strumming,  result 
in  increased  fatigue,  larger  hydrodynamic  forces  and  amplified  acoustic  flow  noise,  and  sometimes  lead 
to  structural  damage  and  eventually  to  costly  failures.  Flow-excited  oscillations  very  often  are  a  critical 
factor  in  the  reliable  and  economical  design  of  underwater  cable  arrays,  mooring  systems,  drilling  risers, 
and  offshore  platforms.  Many  components  of  these  complex  structures  typically  have  a  cylindrical 


cross-section  which  is  conducive  to  vortex  shedding  as  water  flows  past  them. 

As  the  state  of  the  ocean  engineering  art  steadily  progresses,  more  and  more  stringent  demands 
are  being  placed  upon  the  performance  of  cable  structures  and  moorings.  In  particular,  displacement 
tolerances  and  constraints  in  response  to  currents  are  becoming  more  stringent;  fatigue  is  becoming  an 
important  design  consideration;  and  the  sensitivity  of  acoustic  sensors  has  increased  to  the  point  that 
these  sensors  cannot  differentiate  between  legitimate  acoustic  targets  and  slight  variations  in  their  verti¬ 
cal  position.  All  of  these  are  problems  that  are  aggravated  by  cable  strumming.  In  order  to  design  a 
structure  or  cable  system  to  meet  the  constraints  imposed  by  operational  and  environmental  require¬ 
ments.  the  engineer  must  be  able  to  assess  the  effect  of  strumming  on  the  structure  in  question 


1.3  Scope  of  the  Report 

One  purpose  of  this  report  is  to  describe  the  field  test  program  and  to  summarize  me  results  from 
it.  Time  histories  of  the  measured  hydrodynamic  drag  coefficients,  current  speeds,  and  cable  strum¬ 
ming  responses  are  presented  and  discussed  Predictions  are  made  of  the  hydrodynamic  drag  on  a  bare 
cable  and  these  predictions  arc  compared  with  the  field  test  data  for  selected  conditions  when  the  cable 
was  observed  to  be  resonantly  strumming.  The  natural  frequencies  for  the  cable  with  .miched  masses 
were  measured  in  air  and  in  water.  Mode  shapes  were  measured  in  air  Data  for  me  test  runs  con¬ 
ducted  at  the  field  sue  are  presented  in  this  report. 

Validation  of  the  NATFREQ  code  is  an  objective  of  this  report  Calculations  of  cable  natural  fre¬ 
quencies  and  mode  shapes  using  the  NATIREQ  code  were  made  at  NRL  fir  all  of  the  field  test  runs 
A  comparison  is  made  in  this  report  between  the  calculated  natural  frequencies  and  mode  chapes  and 
those  experimental  test  runs  that  have  been  analyzed  in  sufficient  detail  Also,  recommendations  are 
made  of  modifications  that  will  enhance  the  utility  of  and  case  of  access  to  the  code  lor  prospective 
users 
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2.  THE  NATFREQ  COMPUTER  CODE 


2.1  Background  of  NATFREQ 

The  NATFREQ  code  was  developed  ai  lhc  California  Insiiiuic  of  Technology  for  (he  Naval  Civil 
Engineering  Laboratory  (NCEL)  to  calculate  the  natura  frequencies  and  mode  shapes  of  taut  cables 
with  large  numbers  of  attached  discrete  masses  (4).  An  j.gorithm  was  developed  to  solve  the  transcen¬ 
dental  equations  of  cable  motion  by  an  iterative  technique  which  permits  the  calculation  of  extremely 
high  mode  numbers.  The  algorithm  has  been  intpleme  v.ed  as  a  FORTRAN  computer  code  that  is 
available  from  NCEL  and  NRL.  The  algorithm  that  forms  the  basis  of  the  NATFREQ  code  is 
described  in  this  section  of  the  report  A  listing  of  the  code  is  given  in  Appendix  A.  The  methodology 
embodied  in  the  NATFREQ  code  is  computationally  efficient  and  shows  excellent  convergence  proper¬ 
ties  even  for  high  mode  numbers  The  flexibility  of  the  code  provides  the  capacity  for  treating  a  wide 
range  of  cable  system  configurations 

The  static  and  dynamic  analyses  of  a  cable  system  that  experiences  environmental  loading  require 
the  calculation  of  the  hydrodynamic  drag  forces  Strumming  of  the  cables  due  to  vortex  shedding 
increases  the  overall  mean  drag  force  and  causes  a  corresponding  increase  in  the  cable  drag  coefficient. 
The  methods  that  now  have  been  developed  to  calculate  the  drag  due  to  strumming  require  as  inputs 
the  natural  frequencies  and  mode  shapes  of  the  cable  The  NATFREQ  code  provides  a  fast  and  accu¬ 
rate  means  for  computing  these  parameters  for  taut  cables  with  complex  distributions  of  attached 
masses.  It  now  is  possible  with  NATFREQ  to  calculate  the  strumming  drag  on  the  cable  according  to 
the  method  developed  at  NRL  (see  Ref.  1),  including  the  drag  on  the  attached  masses  (sec  Appendix 
C).  These  drag  amplification  factors  can  be  used  as  inputs  to  the  DESADE  and  DECEL1  cable  struc¬ 
ture  analysis  codes  which  are  described  in  Ref  1. 

2.2  Strumming  of  Multi-Segment  Cable  Systems 

The  natural  frequencies  and  mode  shapes  of  cable  oscillation  are  obtained  by  an  iterative  substitu¬ 
tion  algorithm  which  finds  the  solution  satisfying  the  mt  posed  boundary  conditions  of  the  problem 
The  following  development  of  the  solution  is  essentially  the  same  as  the  previous  work  of  Sergev  and 

4 


Iwan  (4) 


Cable  Dynamics.  The  system  considered  is  shown  in  Fig.  1.  It  consists  of  n  cable  segments 
attached  to  n  -  1  masses  The  cable  segments  have  an  effective  mass  per  unit  length  m,  and  a  tension 
T,  which  is  assumed  to  be  constant  over  the  length  /,  of  the  segment.  The  effective  attached  mass 
(including  added  mass)  is  M,. 


Equar&ns  of  Motion.  The  equation  of  motion  for  the  displacement  of  the  Ah  cable  segment, 
assuming  no  bending  rigidity,  is  (see  Fig.  2) 


m, 


92y, 


d2 


y, 


dt1  '  dx2 ' 
The  harmonic  solution  of  this  equation  has  the  form 


(1) 


>,<*,/)  =  Y,(x)eJ’“'  *  (2) 

where  Yfx)  gives  the  shape  of  the  deformation  and  ai  is  the  frequency  of  oscillation.  Substituting  Eq. 

(2)  into  Eq.  (1)  and  solving  the  resulting  equation  gives 

Y,ix I  =  A,  sin  a,iox  4-  B,  cos  a,a>.r.  <  x  <  /, 

t  <31 

=  1 . n 

where 


a,  -  yfmJT,.  <41 

Boundary  Conditions  In  addition  to  satisfying  the  equation  >f  motion,  the  deflection  of  each  cable 
segment  must  satisfy  certain  boundary  conditions.  These  resui  mint  the  geometric  conditions  imposed 
on  the  end?  of  the  cable  assembly,  the  continuity  of  displace  me  at  each  attached  mass,  and  the  bal¬ 
ance  of  forces  at  each  attached  mass. 


Continuity  of  Displacement.  The  displacement  must  be  o  .  nuous  at  each  attached  mass  There¬ 
fore. 

Yflf  <5) 

Substituting  from  Eq.  <31  yields 

B,+ 1  ”  A,  sin  a, u»  /,  +  B,  cos  -  w  I, .  (6) 

If  A,  and  3  are  known.  Eq.  (61  may  be  used  to  find  Bl+]. 
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Force  Bokrnce.  Figure  3  shows  the  forces  acting  at  each  attached  mass.  The  balance  of  forces  at 


each  of  these  points  gives  for  small  angles 


T,.i 


Byl+t 

Bx 


-  r 


By, 

Bx  [»-/, 


M, 


B2yl+ 1 

dr 2 


U-o 


Substituting  from  Eqs  (2)  and  (3)  yields 


(7) 


7^,0  -  'a, A,  cos  a, to  /, 

—  a,B,  sin  a,(u  /, )  T,  =  —  M,<u  B , + 1 

/  -  I . n. 

Using  Eq.  (6)  and  solving  Eq.  (8)  for  A,+1  gives 


(8) 


,4,. [  =  — — - —  (fa,  T,  cos  a,w  /,  —  Midi  sin  a,oi  I,)  A 

-  la,  7",  sin  a,di  I.  +  M,di  COS  a,co /,)«,].  (9) 

If  A,  and  B  are  known.  Eq.  (9)  may  be  used  to  find  A,  + 1 

Geometric  Boundary  Conditions.  At  the  left  hand  end  ol  the  cable  assembly,  the  displacement  is 
assumed  to  be  zero.  Thus. 

K,(0)  =  0.  <101 

Substituting  from  Eq  13)  thus  implies  that 


0,  =  0.  Ill* 

Since  the  scale  of  the  deflected  shape  of  (he  cable  (mode  shape)  is.  at  this  point,  arbitrary,  let 


A,-  I.  (12) 

With  conducts  <  1 1  •  and  <121  on  1 1  and  B j  and  Eqs.  (6)  and  (9)  lor  4,*|  a  id  B,.  |.  all  subsequent  A  s 

and  B' s  can  be  determined  provided  m  (the  natural  frequency)  is  known 


The  system  must  satisfy  one  additional  boundary  condition  at  the  right  hand  end  of  the  cable 
assembly  ,  w  ere  the  displacement  is  again  assumed  to  be  /ero.  This  gives 

Y„U,)  =  0-  K,irl«!>.  (13) 

Equation  1 1 ' '  in  turn  implies  that 

«„  m-0.  1141 
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The  values  of  <u  which  give  solutions  satisfying  condition  (14)  are  the  natural  frequencies  of  the  sys¬ 
tem  This  is  a  shooting  technique. 


2.3  Solution  Algorithm  for  Mode  Shapes  and  Frequencies 

If  at  is  varied  from  zero  to  some  large  value  and  the  corresponding  values  of  S„+|  calculated,  the 
result  will  be  as  shown  in  Fig  4.  F.ach  point  for  which  B„  +  l  —  0  represents  a  valid  solution  of  the  free 
oscillation  problem.  The  ca*  so  obtained  are  the  natural  frequencies  of  the  system.  Mathematically 
•h?re  are  an  infinite  number  of  such  frequencies. 

The  mode  shape  associated  with  each  natural  frequency  <ok  will  be  denoted  by  F,u)(ar).  Then, 

=  /1,<U  ski  a,oikx  -I-  B,ik'  cos  a,wkx.  (15) 

Let  <jj.  be  the  Aih  natural  frequency.  Then  the  deflected  shape  of  the  cable  system  will  be  such  that  the 
number  of  internal  zero  crossings  (nodes)  is  equal  to  k  -  1.  The  mode  number  of  a  particular  mode 
srape  may  therefore  be  determined  by  counting  the  number  of  internal  zeros  associated  with  the  func- 
'uin  Y'k> .  (  —  l . ti. 

One  of  the  example  cases  computed  with  NATFREQ  was  a  15400  ft  long  cable  with  380  attached 
bodies  The  calculated  modal  pattern  for  mode  number  162  is  shown  in  Fig.  5.  The  complexity  of  the 
response  is  evident.  Some  additional  examples  and  a  comparison  of  NATFREQ  results  with  a  small- 
scaie  laboratory  simulation  are  given  by  Sergev  and  Iwan. 

2.4  Summary  of  the  Solution  Procedure 

The  solution  process  for  the  mode  shapes  and  frequencies  can  be  summarized  as  follows  (4): 

1  Assume  a  value  for  oj. 

2  I  ci  fl,  -  0.  .4,  *  I 
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3.  Solve  for  fl2,  A2\  B},  A 3,  . ..  ;  B„,  A„:  from  Eqs.  (6)  and  (9). 


4  Check  for  Bn> |  —  0.  If  #„tl  ^  0,  compare  .<nh  previous  value  and  estimate  a  new  trial 
value  for  a>. 

5  Go  to  step  2  and  repeat  until  Bnt[  is  less  than  some  prescribed  value  or  the  change  in  <«t 
is  less  than  some  prescribed  limit 

6  Determine  the  mode  number  by  calculating  the  number  of  internal  zeros  of  the  mode 
shape  K,ul(.v);  i  -  1.  . ..  n 

A  listing  of  the  NATFRF.Q  computer  code  presently  in  use  at  NRL  is  given  in  Appendix  A. 
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3.  THE  TEST  SITE  AND  INSTRL  MENTATION 


3.1  The  Test  Site 

The  site  chosen  for  the  experiment  was  a  sandbar  located  at  the  mouth  of  Holbrook  Cove  near 
Castine,  Maine.  This  was  the  same  site  used  for  previous  experiments  during  the  1970's  by  Vandiver 
and  Mazel  (5.6).  At  low  tide  the  sandbar  was  exposed  allowing  easy  access  to  the  test  equipment, 
while  at  high  tide  it  was  covered  by  about  ten  feel  of  water.  The  test  section  was  oriented  normal  to 
the  direction  of  the  current  which  varied  from  0  to  2.4  ft/s  over  the  tidal  cycle  with  only  small  spatial 
variation  over  the  test  section  length  at  any  given  moment. 

The  data  acquisition  station  for  the  experiment  was  the  R/V  Edgerton  which  was  chartered  from 
the  MIT  Sea  Grant  Program  The  Edgerton  was  moored  for  the  duration  of  the  experiment  approxi¬ 
mately  300  feet  from  the  sandbar  and  connected  to  the  instruments  on  the  sandbar  by  umbilicals. 

Prior  to  the  data  acquisition  phase  of  the  experiment,  several  days  were  needed  to  prepare  the 
site  A  foundation  for  the  experiment  was  needed  to  anchor  the  supports  which  were  to  hold  the  ends 
of  the  test  cylinders.  To  accomplish  this,  six  4.5-inch  diameter  steel  pipes  were  water  jetted  into  the 
sandbar  utilizing  the  fire  pump  aboard  the  Edgerion.  These  six  pipes  were  made  of  two  five  foot  sec¬ 
tions  joined  by  couplings  so  that  the  overall  length  of  each  was  ten  feet.  In  addition,  one  two-inch 
diameter  by  six-foot  long  steel  pipe  was  jetted  into  the  sandbar  to  be  used  as  a  current  meter  mount. 
Finally,  a  section  of  angle  iron  was  clamped  to  the  pipe  used  to  support  the  drag  measuring  mechanism 
and  attached  to  another  support  pipe  to  prevent  any  rotation  of  the  drag  mechanism  mount.  Figure  6 
shows  a  schematic  diagram  of  the  set-up  of  the  experiment. 

3.2  Test  Instrumentation 

Drag  measurement  system  The  drag  measurement  system  was  located  at  the  west  end  of  the  cable 
system  as  shown  in  Fig  6  The  device  was  welded  onto  a  support  pipe  2.5  feet  above  the  mud  line 
The  mean  drag  force  at  the  termination  of  the  cable  was  used  to  generate  a  moment  about  a  Irecly 


MODE  NO  •  162 


l-ig  5  —  The  NATFRF.Q-computed  mode  shape  (mode-number  n  —  162 )  for  a  15400  fl  long  marine  cable  wiih  380  attached 
masses;  from  Sergcv  and  Iwan  (4)  1  in  “  25.4  mm. 


routine  vertical  shaft  located  a  few  inches  beyond  the  termination  point.  The  bearings  supporting  the 
shaft  carried  the  entire  tension  load  without  preventing  roution.  The  moment  was  balanced  by  a  load 
cell  which  restrained  a  lever  arm  connected  to  the  shaft  (see  Fig.  7).  From  the  known  lever-arm 
lengths  and  the  load  cell  measurements  the  mean  drag  force  on  one  half  of  the  cable  could  be  deter¬ 
mined.  The  lead  cell  signal  was  carried  by  wires  in  the  cable  and  umbilical  to  the  Edgerton  where  it  was 
conditioned  and  recorded. 

Current  measurement  system.  The  current  was  measured  by  a  Neil  Brown  Instrument  Systems 
URCM-2  Acoustic  Current  Meter  located  12.5  ft  from  the  west  end  of  the  test  cable  and  2  ft  upstream. 
Ii  was  set  so  that  it  determined  both  normal  and  tangential  components  of  the  current  at  the  level  of 
th :  test  cable  Signals  from  the  current  meter  traveled  through  umbilicals  to  the  Edgerton  where  they 
were  monitored  and  recorded.  In  addition,  a  current  meter  traverse  was  made  using  an  Endeco  current 
meter  to  determine  any  spatial  variations  in  current  along  the  test  section  The  current  was  found  to  be 
spatially  uniform  to  within  ±3.0  percent  from  end  to  end  for  all  but  the  lowest  current  speeds  (V  < 
0  5  ft/s). 

Tension  measurement  system  The  tension  measuring  and  adjusting  system  was  located  at  the  east 
end  of  the  experimental  test  set  up  (see  Fig.  6).  Extensions  were  maJe  to  the  two  inner  water-jetted 
posts  at  this  end.  As  shown  in  the  diagram,  a  five-foot  extension  was  made  to  the  center  post  and  a 
three  foot  extension  was  made  to  the  innermost  post.  This  three-foot  extension  was  different  from  the 
rest  in  that  its  attachment  to  the  jetted  pipe  at  the  mudline  was  a  pin  connection  as  compared  to  the 
standard  pipe  couplings  used  on  the  other  extensions.  A  hydraulic  cyl  nder  was  mounted  2.5  ft  above 
the  mudline  onto  this  pivoting  post.  The  lest  cable  in  the  experiments  was  connected  at  one  end  to  this 
hydraulic  cylinder  and  at  the  other  end  to  the  drag  measuring  device.  To  the  back  of  the  hydraulic 
cylinder  one  end  of  a  Sensotec  Model  RM  In-Line  load  cell  war  connected.  The  other  end  of  the  cell 
was  attached  s  la  a  cable  to  the  center  post.  The  output  from  the  tension  load  cell  was  transmitted 
through  the  umbilicals  lo  the  Edgerton  where  it  was  monitored  Hydraulic  hose  ran  from  a  pump  on  the 
Edgerton  to  the  hydraulic  cylinder  so  that  the  tension  could  be  changed  as  desired  Additional  details 
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Fig  8  —  A  cross-section  of  the  test  cable,  from  Vandiver  and  Griffin  (7). 
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concerning  the  test  instrumentation  are  given  by  McGlothiin  (8). 


3.3  Data  Acquisition  Systems 

During  the  experiment,  data  taken  from  the  instruments  on  the  sandbar  were  recorded  in  two 
ways.  First,  analog  signals  from  the  fourteen  accelerometers  in  the  cable  as  well  as  current  and  drag 
were  digitized,  at  30.0  Hz  per  channel,  onto  floppy  disks  using  a  Digital  Equipment  MINC-23  Com¬ 
puter.  Second,  analog  signals  from  the  drag  cell,  current  meter,  and  six  accelerometers  were  recorded 
by  a  Hewlett-Packard  Model  3968A  Recorder  onto  eight-track  tape.  The  disks  were  limited  to  record 
lengths  of  eight  and  one  half  minutes  and  were  used  to  take  data  several  times  during  each  two  and  one 
half  hour  data  acquisition  period.  A  Hewlett-Packard  3582 A  Spectrum  Analyzer  was  set  up  to  monitor 
the  real  time  outputs  of  the  accelerometers.  The  eight-track  tape  was  used  to  provide  a  continuous 
record  of  the  complete  two  and  one  half  hour  test  run 

4.  THE  TEST  CABLE  SYSTEM 

4.1  The  Cable 

A  75-foot  long  composite  cable  was  developed  at  MIT  specifically  for  the  experiments  that  were 
conducted  in  the  summer  of  1981.  Figure  8  shows  a  crrss-section  of  the  test  cable.  The  outer  sheath 
for  this  cable  was  a  75-foot  long  piece  of  clear  flexible  PVC  tubing,  which  wa>  1. 25-in.  O.D  by  1  0-in 
I  D.  Three  0. 1 25-in.  stainless  steel  cables  ran  through  the  tubing  and  served  as  the  tension  carrying 
members.  A  cylindrical  piece  of  0.5-in.  O.D.  neoprene  rubber  was  used  to  keep  the  stainless  steel 
cables  spaced  120  degrees  apart.  The  neoprene  rubber  sr.jcer  was  continuous  along  the  length  except  at 
seven  positions  where  biaxial  pairs  of  accelerometers  were  placed.  Starting  at  the  east  end.  these  posi¬ 
tions  were  at  L/8,  L/6,  L/4.  2L/5.  L/2,  SL/8,  and  3L.  -  These  accelerometers  were  used  to  measure 
the  response  of  the  cable  as  it  was  excited  by  the  vortex  shedding  The  accelerometers  were 
Sundstrand  Mini-Pal  Model  2180  Servo  Accelerometers  wr.ich  were  sensitive  to  the  direction  of  gravity 
The  biaxial  pairing  of  these  accelerometers  made  it  r  'ssible  to  determine  their  orientation  and  to 


extract  real  vertical  and  horizontal  accelerations  of  the  uMc  at  the  seven  locations 


Three  bundles  of  ten  wires  each  ran  along  the  sides  of  the  neoprene  spacer  to  provide  power  and 
signal  connections  to  the  accelerometers  and  also  to  provide  power  and  signal  connections  to  the  drag 
measuring  system.  Finally,  an  Emerson  and  Cuming  flexible  epoxy  was  used  to  fill  the  voids  in  the 
cable  and  make  it  watertight.  The  weight  per  unit  length  of  this  composite  cable  was  0  77- lb/ft  in  air. 
A  photograph  of  the  test  cable  being  installed  for  one  of  the  test  runs  is  shown  in  Fig  9. 

4.2  The  Attached  Masses 

In  some  experiments,  lumped  masses  were  fastened  to  the  bare  cable  to  simulate  the  effects  of 
sensor  housings  and  other  attached  bodies.  The  lumped  masses  were  made  of  cylindrical  PVC  stock 
and  each  was  12.0-in.  long  and  of  3.5-in.  diameter.  A  1.25-in.  hole  was  drilled  through  the  center  of 
each  lumped  mass  so  that  the  cable  could  pass  through.  In  addition,  four  0.625-in.  holes  were  drilled 
symmetrically  around  this  1.25-in.  center  hole  so  that  copper  tubes  filled  with  lead  could  be  inserted  to 
change  the  mass  of  lumps.  In  the  field,  it  was  difficult  to  force  the  cable  through  the  holes  drilled  in 
the  PVC  so  the  masses  were  split  in  half  along  the  length  of  their  axes.  The  masses  were  then  placed 
on  the  cable  in  halves  and  held  together  by  hose  clamps.  Different  tests  were  run  by  varying  the 
number  and  location  of  lumped  masses  and  by  changing  the  mass  of  the  attachments.  A  photograph  of 
one  of  the  masses  attached  to  the  test  cable  is  shown  in  Fig.  10. 

5.  THE  TEST  MATRIX  AND  SPECIFICATIONS 

5.1  Test  Matrix  and  Specifications 

The  test  cases  conducted  for  the  combination  of  the  cable  and  Ihe  various  attached  masses  is 
given  in  Table  I  The  input  data  for  the  NATFREQ  calculations  of  the  twenty  test  runs  are  listed  in 
Appendix  B.  The  specifications  for  the  cable  and  the  various  attached  mass  configurations  are  given  in 
Table  2.  The  notation  for  the  location  of  and  type  of  attached  mass  in  Table  1  is.  as  one  example. 

T 

This  notation  means  that  the  mass  was  heavy  (//),  i  c..  that  lead  weights  were  inserted  into  the  PVC 
cylinder,  and  that  the  attachment  point  was  located  5/8  of  the  cable  length  from  the  reference 


17 


Table  1  —  Test  Run  Sequence 
1981  Castine  Bay  Field  Experiments, 
Marine  Cables  with  Attached  Masses 


Run 

Location,  Type  of 

Attached  Mass 

Tension 

(pounds) 

Medium 

1 

j  (L),  ( L )» 

512— 

Airf 

2 

same  as  above 

450 

Water 

3 

~  (L).  y  (L),  —  (L),  ~  (L),  ~  (I)  ^ 

500 

Air 

4 

same  as  above 

520 

Water 

5 

j  (L).  ^  ( L ),  ( L ),  ^  U) 

532 

Air 

6 

same  as  above 

475 

Water 

7 

4  (L).  4  ( L ).  ^  (L),  4  (L).  <£)•  <i>-  HT  <*> 

8  4  8  2  8  4  o 

504 

Air 

8 

same  as  above 

475 

Water 

9 

j  (H)  ^  (H) 

540 

Air 

10 

same  as  above 

500 

Water 

11 

H >.  |-( H ),  y(//),  -y(//),  ^4  (tf) 

728 

Air 

12 

same  as  above 

580 

Water 

13 

j(H),  ^(//),  -y^(//).  ^  ( H ) 

650 

Air 

14 

same  as  above 

650 

Water 

15 

y(tf).  y<//).  ^  (H) 

732 

Air 

16 

same  as  above 

556 

Water 

17 

4  (H).  4  W).  ~4  (#>,  4  -T1  <£> 

6  2  6  3 

800 

Air 

18 

same  as  above 

765 

Water 

19 

4  <L>.  4  (D-  4  <#/).  (//).  ~  (//•■  4r 

6  2  3  8  4  8 

800 

Air 

20 

same  as  above 

800 

Water 

Ml  tensions  were  uken  from  data  records  from  which  natural  frequencies  were  obtained. 
'</.)  "*  light  mass.  < H)  —  Heavy  mass 
“One  pound  -  4  4  Newtons  (SI  units) 
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Table  2  —  Cable  and  Attached  Mass  Parameters 
1981  Castine  Bay  Field  Experiments, 
Marine  Cables  with  Attached  Masses 


Cable  Specifications: 

L  -  75.0  feet 

D  “  1.25  inches  in  diameter 
w  =*  0.7704  pounds/foot  in  air 
n  -  0.0239  slugs/foot  in  air 

SGC  *■  1 .408  specific  gravity 

i  ' 

PVC  Attached  Mass  Specifications: 

Length  —  12  0  inches 

Outside  diameter  —35  inches 

Inside  diameter  —  1  25  inches 

Ballast  hole  diameter  «  0  625  inches 

Weight  in  air  — 

i  ' 

1  ' 

No  lead  ballast,  4.4  pounds  (2  kg) 

Weight  of  water  trapped  in  the  5/8  inch  ballast 
holes  when  submerged  in  water,  0  55  pounds  (0.25  kg) 
Weight  with  lead,  9.97  pounds  (4.52  kg) 

t 

w 


lermmation  point  The  notation  L  in  parentheses  (£)  means  that  the  light  mass  (no  lead  weights)  was 
attached  at  that  location.  The  two  specific  gravities  were  chosen  to  match  the  cable's  specific  gravity  SG 
=  13  to  14  (nominal)  and  the  specific  gravity  SG  =  2  (nominal)  of  typical  Navy  hydrophones  and 
other  cable-mounted  sensor  housings. 

5.2  Slaek  Cable  Limitations 

Extensible  slack  cables  are  characterized  by  complex  dynamic  response  behavior  that  is  dependent 
upon  the  sag-to-span  ratio  and  the  elastic  properties  of  the  cable.  One  especially  complicated  feature  of 
ti'is  response  is  a  frequency  "crossover*  phenomenon.  This  modal  crossover  is  a  complex  phenomenon 
associated  with  the  dynamics  of  slack  cables  with  small  sag-to-span  ratios.  At  the  crossover  three  modes 
of  the  cable  have  the  same  natural  frequency  and  include  a  symmetric  in-plane  mode,  an  anti¬ 
symmetric  in-plane  mode  and  an  out-of-plane  or  sway  mode.  The  symmetric  modes  contain  an  even 
number  of  nodal  points  along  the  cable  while  the  anti-symmetric  modes  contain  an  odd  number  of 
nodes.  The  dynamics  of  slack  marine  cables  are  discussed  further  in  Appendix  B  of  Ref.  1,  and  in 
several  other  references  cited  there. 

The  onset  of  slack  cable  effects  occurs  near  the  limit  of  cable  sag-to-span  s/J  — -  0.  and  the  critical 
tension  Han  can  be  estimated  with  the  equation  (see  Ref  1 ) 

Han  =  0.93  iWlExA,)Vi.  (16) 

where 


W  —  total  cable  weight; 

E.  =  cable  elastic  modulus; 

* 

A ,  *  cable  cross-section  area. 

The  weight  of  the  cable  in  water  is 

W ■-  (p.eJL.^dSG.  -  1)  (17) 

w  here 
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w  All 


'TtJr*  .«* 


SGC  —  cable  specific  gravity; 

Lc  —  cable  length; 

Pw  8  water  density,  gravitational 
acceleration 

The  Castine  Bay  test  cable  (Lc  “  75-ft,  Ac  =*  1,23-in.2)  in  sea  water  has  the  following  total  net  weight 
(in  water), 

W  -  41.0  (SGC  -  1)  lb. 

For  a  cable  with  a  nominal  value  of  specific  gravity  SGf  —  1.41.  W  —  16.8  lb.  The  product  ECAC 
was  measured  at  MIT  for  the  Castine  test  cable.  The  result  was  £r.4..  —  2(105)  lb.  Based  upon  this 
value  the  critical  tension  in  water  is 

//cm  =  0.93  (ECAC  B/2)i/-'  lb  -  360  lb. 

From  this  estimate  a  minimum  tension  H  =  400  lb  was  recommended  for  the  Castine  Bay  field  experi¬ 
ments  in  order  to  minimize  slack  cable  effects.  The  NATFREQ  code  algorithm  is  valid  only  for  taut 
cables  v*  ith  arrays  of  attached  masses. 

6.  NATFREQ  PREDICTIONS 

The  NATFREQ  code  was  used  to  simulate  all  of  the  test  runs  listed  in  Table  1.  All  of  the  com¬ 
putations  were  done  on  the  Hewlett-Packard  Model  1000  computer  in  the  Fluid  Dynamics  Branch  at 
NRL.  A  listing  of  the  code  is  given  in  Appendix  A.  The  listed  version  of  the  code  is  identical  to  that 
discussed  by  Sergev  and  Iwan  (4)  except  for  some  minor  input/output  modifications  Input  to  the  code 
typically  consists  of  cable  properties  (density,  tension,  drag  coefficient,  etc  ),  attached  mass  properties 
(density,  cross-section  area,  drag  coefficient,  etc.),  and  cable  segment  length  (point  masses  are 
assumed)  as  shown  in  Appendix  B.  The  added  mass  of  both  the  cable  and  the  attached  masses  was 
accounted  for  in  all  of  the  in-water  test  simulations.  Some  elaborations  on  the  algorithm  employed  in 
NATFREQ.  but  not  included  in  the  code,  are  given  by  Iwan  (9). 

A  summary  of  the  computer  simulation  is  given  here.  Typically  only  the  first  six  or  seven  cable 


modes  were  excited  by  the  tidal  current  regime  at  the  Castine  field  site  However,  in  some  cases  cable 
modes  up  to  n  **  11  have  been  identified.  The  computed  natural  frequencies  and  mode  shapes  for  five 
pairs  of  typical  test  runs  are  plotted  in  Figs.  11  to  18.  The  configurations  chosen  for  the  plots  are 
representative  of  the  field  test  measurements  that  are  discussed  in  the  next  section.  Three  cases  with 
light  masses  attached  (PVC  only)  and  two  cases  with  heavy  masses  (PVC  plus  lead  weight  inserts)  are 
plotted  in  the  figures 

Several  points  concerning  the  natural  frequencies  are  of  note.  First,  the  added  mass  effect  of  the 
water  is  clearly  evident  in  all  of  the  computed  frequencies.  As  expected,  the  added  mass  due  to  the 
cable  vibration  reduces  the  frequency  of  any  given  in-water  mode  below  the  corresponding  in-air  modal 
frequency.  An  added  mass  coefficient  of  Cam  “  1  was  assumed;  i.e.,  the  added  masses  of  the  cable  and 
cylindrical  attachments  are  equal  to  the  respective  displaced  masses  of  water.  This  point  is  discussed 
further  in  a  following  section  of  this  report. 

Second,  the  distribution  of  computed  modal  frequencies  is  dependent  upon  the  number  and  spac¬ 
ing  of  attached  masses.  For  evenly  spaced  masses  there  are  discontinuous  jumps  in  the  natural  fre¬ 
quencies  at  multiples  of  the  number  of  cable  segments  This  is  clearly  evident  in  Fig.  1 1  (two  and  five 
attached  masses)  and  the  same  effect  can  be  observed  in  Figs.  14,  16  and  18  (five,  six  and  seven 
attached  masses).  The  discontinuity  emerges  when  the  mode  number  n  is  equal  to  the  number  of  cable 
segments  or  a  multiple  thereof.  The  difference  in  the  cable  vibration  pattern  due  to  different  numbers 
of  attached  masses  is  evident  from  the  N  ”  2  and  N  “  5  results  in  Fig  11  The  discontinuities  occur 
at  n  “  3,  6  and  9  for  the  cable  with  two  masses,  whereas  the  first  discontinuity  occurs  at  n  «  6  for  the 
cable  with  five  masses.  Note  that  the  n  **  6  natural  frequencies  (and  mode  shapes)  are  equal  for  the 
two  cases.  This  behavior  also  was  computed  by  Sergev  and  Iwan  and  verified  by  their  small  scale  exper¬ 
iment  with  a  thin  wire  and  four  or  six  attached  masses  in  air  (4). 

There  is  a  relatively  simple  physical  explanation  for  the  discontinuous  behavior  of  the  natural  fre¬ 
quencies.  Consider  as  an  example  the  computed  frequencies  in  Fig.  14,  for  test  runs  seven  and  eight 
with  seven  light  masses.  For  the  first  seven  mode  shapes  plotted  in  Figs.  1 S (a)  and  15(b)  some  or  all 


of  the  masses  are  in  motion  at  their  various  locations  in  the  vibration  pattern.  Thus  the  natural  fre¬ 
quencies  of  the  cable  are  reduced  from  the  bare  cable  values  due  to  the  inertia  of  the  masses.  When 
the  vibration  frequency  reaches  the  n  —  8  mode  all  of  the  seven  evenly-spaced  attached  masses  are 
located  at  nodes  of  the  modal  vibration  pattern.  Then  the  mode  shape  and  frequency  are  identical  to 
those  of  the  bare  cable.  The  natural  frequencies  for  the  n  —  9  and  10  modes  fall  below  the  bare  cable 
frequencies  when  the  masses  attached  to  the  cable  again  are  in  motion.  Similar  behavior  can  be 
observed  in  all  of  the  natural  frequencies  plotted  in  Figs.  11,  14,  16  and  18. 

The  distribution  of  six  masses  for  the  test  runs  plotted  in  Fig.  16  is  unevenly  spaced.  The  discon¬ 
tinuity  between  the  n  —  6  and  n  -  7  modes  still  is  clear,  but  the  n  -  7  mode  in  air  frequency  is  not 
equal  to  the  bare  cable  value  as  in  the  runs  discussed  previously.  The  mode  shapes  plotted  in  Fig.  17 
show  that  several  of  the  masses  are  in  motion  for  n  -  7,  so  that  both  the  in-air  and  in-water  natural 
frequencies  are  less  than  the  bare  cable  value.  There  is  a  difference  in  tension  between  the  two  tests 
but  that  has  a  minimal  effect  on  the  major  features  of  the  overall  cable  system  response. 

Two  otherwise  identical  runs  with  five  light  and  five  heavy  masses  are  plotted  in  Fig.  18.  Test 
runs  11  and  12  were  recomputed  with  the  cable  tension  T  —  500  pounds  to  match  the  tensions  of  runs 
3  and  4.  The  inertial  effect  of  the  heavier  of  the  two  attached  masses  is  evident  from  the  reduced 
natural  frequencies  in  air  and  in  water  below  the  corresponding  values  for  the  lighter  masses.  The 
masses  are  evenly  spaced  and  the  discontinuity  between  the  n  -  5  and  n  -  6  modes  appears  as 
expected.  For  the  n  -  6  mode  both  sets  of  masses  are  at  the  same  nodal  points  of  the  vibration  pattern 
and  they  have  no  effect.  The  two  respective  pairs  of  natural  frequencies  in  air  and  in  water  are  equal. 


NATURAL  FREQUENCY,  fn/Hz 


Fig.  11  —  ln-«ir  ind  in-water  natural  frequencies  for  a  cable  with  attached  masses,  computed  with  NATFREQ.  — .  Test  Run  1. 
t.i.u  1  — .  Test  Run  2; - ,  Test  Run  3; . ,  Test  Run  4;  O,  bare  cable. 
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MODE  NO.  -  2 


MODE  NO.  -  3 


MODE  NO.  -  4 


Fig  12(a)  —  Test  Run  I.  Table  1 

Fig.  12  —  Mode  shapes  for  a  cable  with  two  light  attached  masses,  computed  with  NATFREQ. 
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1981  CASTINE  BAY  TESTS  *  MARINE  CABLES 
WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  2 


Fig  11(b)  —  Continued 


0 


31 
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MODE  NO. 
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1981  CASTINE  BAY  TESTS  •  MARINE  CABLES 
WITH  ATT  ACHED  SENSOR  HOUSINGS  •  RUN  3 


Ftg  13(«l  —  Continued 
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1981  CASTINE  BAY  TESTS  *  MARINE  CABLES 
WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  4 


Fig  13(bl  —  Conunucd 


ATTACHEO  MASS  LOCATIONS: 
L/8  _  3L/8  5L/8  7L/8 

LA  L/2  3L/4 


m-water  natural  frequencies  for  a  cable  with  seven  light  attached  masses,  computed  with  NATFRF.Q 


-0.0 


7S.  0 


NATURAL  FREQUENCY,  fn/Hz 


Fig.  16  —  In-air  and  in-water  natural  frequencies  for  a  cable  with  six  heavy  attached  misses,  computed  with  NATFREQ  — . 
Run  15  Table  1;  — —  Run  16;  0.  bare  cable 


40 


1981  CASTINE  BAY  TESTS  *  MARINE  CABLES 
WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  15 

a.  a  _  MODE  NO.  -  1 


Fig  ITtal  —  Test  Run  15.  Table  I 

Ft*  17  —  Mode  shapes  for  a  cable  with  six  heavy  attached  maaaev  computed  with  NATFREQ 


1981  CASTINE  BAY  TESTS  »  MARINE  CABLES 
WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  16 


a.  • 


MODE  NO.  -  5 


1981  CASTINE  BAY  TESTS  *  MARINE  CABLES 
WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  16 


*.  ■ 


MOOE  NO 


9 


NATURAL  FREQUENCY,  fn/Hz 


MODE  NUMBER,  n 

Fig.  18  —  In-air  and  m-water  natural  frequencies  for  a  cable  with  five  attached  masses,  computed  with  NATFREQ.  — ,  Run  3, 
Table  I. - ,  R.n  4. - ,  same  as  Run  11,  except  T  —  500  lb; - ,  same  as  Run  12,  except  T  •  500  lb;  O.  bare  cable 
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The  modal  distributions  of  vibration  amplitude  along  the  cable,  respectively  corresponding  to  the 
plotted  frequencies,  are  shown  in  Figs.  12,  13,  15,  and  17.  The  vertical  scales  are  uniform  but  arbi¬ 
trary,  having  been  chosen  only  for  clarity.  The  maximum  strumming  vibration  amplitude  for  a  given 
cable  mode  is  dependent  upon  the  effective  reduced  damping  parameter  ks\  see  Eq.  (C3)  of  Appendix 
C.  For  vibrations  of  a  cable  in  water  the  typical  level  of  strumming  amplitudes  is  on  the  order  of  (±) 
one  cable  diameter.  Vortex-induced  vibrations  of  cables  and  structures  in  air  are  an  order  of  magnitude 
lower  in  amplitude.  Thus  even  though  the  relative  mode  shapes  in  air  and  in  water  may  appear  to  be 
very  similar  on  a  nondimensional  basis,  the  actual  distributions  of  strumming  amplitude  levels  spanwise 
along  a  cable  are  very  different  in  the  two  medias. 

It  is  evident  from  the  plotted  mode  shapes  that  the  attached  masses  do  not  act  in  general  as  nodes 
of  the  overall  modal  cable  pattern.  This  is  true  for  all  of  the  cases  shown  with  up  iq  seven  attached 
masses.  The  complexity  of  the  modal  vibration  patterns  along  the  cable  increases  in  the  higher  modes. 
Thus  the  higher  mode  shapes  for  a  cable  with  large  numbers  of  attached  masses  cannot  be  guessed 
intuitively  from  previous  experience.  However,  the  mode  shape  and  peak  strumming  amplitude  must 
be  accurately  known  in  order  to  correctly  estimate  the  drag  effects  due  to  strumming  of  a  complex  cable 
segment. 

7.  MEASUREMENTS  OF  CABLE  STRUMMING 
7,1  Natural  Frequencies,  Mode  Shapes  and  Damping 

Measurements  in  Air.  The  natural  frequencies  and  damping  ratios  weie  measured  in  air  for  each 
cable  test  run  with  attached  masses  and  for  several  tests  with  the  bare  cable  (10,  11).  A  method  for 
extracting  the  natural  frequencies,  mode  shapes,  cable  material  damping,  and  strumming  displacement 
from  the  outputs  of  the  accelerometer  pairs  was  developed  and  calibrated  as  part  of  the  data  reduction 
phase  of  the  test  program  (12).  A  summary  of  the  experimental  data  for  in-air  natural  frequencies, 
mode  shapes,  and  damping  ratios  is  given  by  Vandiver  (13).  An  example  of  the  good  agreement  that 
was  achieved  between  the  measured  and  calculated  natural  frequencies  is  given  in  Table  3  The  natural 
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Table  3  —  Natural  Frequencies  and  Damping  of  a  Bare  Cable  (In  Air) 
1981  Castine  Bay  Field  Experiments, 

Marine  Cables  with  Attached  Masses; 

From  Refs.  11  and  13 
(Cable  Tension  -  792  pounds) 


Mode 

Natural  Frequency 

Natural  Frequency 

Material  Damping  Ratio,  ( 

(predicted,  Hz) 

(observed,  Hz) 

Vertical 

based  on 
energy  contained 
in  principal 
direction  (() 

based  on  total 
energy 

1 

1.21 

2 

2.42 

2.40 

0.0065  <  <  0.0084 

0.0024  ^  <  0.0065 

3 

3.63 

3.60 

0.0045  <  0.0049 

0.0017  <  £  <  0.0056 

4 

4.84 

4.67 

U  “  0  005 

0  0021  ^  C  <  0.0046 

5 

6.05 

5.86 

6 

7.26 

7.23 

7 

8.47 

a 

8 

968 

9 

10.89 

10.92 

10 

12.10 

frequencies  were  calculated  with  the  standard  equation 


-  2 ir/,,  - 


(18) 


for  the  undamped  frequencies  of  a  taut  cable  per  unit  length.  Here  m  is  the  virtual  (physical  4-  added) 
mass  of  the  cable  per  unit  length.  Excellent  agreement  was  obtained  for  the  measured  and  calculated 
natural  frequencies  of  the  bare  cable,  and  considerable  confidence  was  gained  in  the  data  reduction 
technique  in  the  process.  The  very  small  material  damping  of  the  cable  is  characteristic  of  marine 
cables  (1 ). 


The  accelerometer  pairs  embedded  in  the  cable  were  of  the  force-balanced  type  and  sensitive  to 
the  direction  of  gravity.  Thus  the  true  horizontal  and  vertical  accelerations  of  the  cable  could  be 
obtained.  Once  the  true  vertical  and  horizontal  accelerations  were  found,  it  was  necessary  to  undertake 
a  complex  process  to  determine  the  displacements.  To  do  this,  a  step-by-step  process  was  developed 

A  summary  of  the  procedure  for  double  integrating  a  digital  acceleration  signal  follows.  The  pro¬ 
cedure  consists  of  the  following  eleven  steps  as  outlined  by  Jong  (12): 

Step  I  Find  the  rotated  angle  of  the  accelerometers  and  perform  vector  rotation  to  recover  real 
vertical  and  horizontal  accelerations. 

Step  2  Fit  a  linear  least-squares  zeio  baseline  to  the  acceleration  signal  to  remove  the  DC  offset 
and  any  linear  trend 

Step  3  Obtain  the  acceleration  spectrum  and  from  it  compute  the  theoretical  nsplacement  spec¬ 
trum,  as  an  aid  in  the  determination  of  the  cutoff  frequency  in  the  high-pass  filter  The 
theoretical  fast  Fourier  transform  (FFT)  of  the  displacement  is  obtained  by  dividing  the 
magnitude  of  the  acceleration  FFT  by  at2 

Step  4  High-pass  filler  the  acceleration  signal  using  an  infinite  impulse  response  (HR)  elliptic  filter 

to  remove  low  frequency  noise  components. 
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Step  5  Integrate  the  acceleration  signal  using  the  Schuessler-Ibler  integrator  to  obtain  the  velocity. 

Step  6  Fit  a  linear  least-squares  zero  baseline  to  the  velocity  signal  to  remove  any  DC  offset  and 

linear  trend. 

Step  7  High-pass  filter  the  velocity  signal  using  an  HR  elliptic  filter  to  remove  any  low  frequency 
noise  components  that  were  expanded  in  step  5. 

Step  8  Integrate  the  velocity  signal  using  the  Schuessler-Ibler  integrator  to  obtain  the  displacement. 

Step  9  Fit  a  linear  least-squares  zero  baseline  to  the  displacement  signal. 

Step  10  High-pass  filter  the  displacement  signal  using  an  IIR  elliptic  filter  to  remove  any  low  fre¬ 
quency  noise  components  that  were  expanded  in  step  8. 


Step  11  Plot  summary  data,  e.g.  root-mean-squares,  spectra,  time  series  of  displacement,  velocity, 
acceleration  and  two  dimensional  cylinder  or  cable  motion  time  series. 

From  the  horizontal  and  vertical  acceleration  spectra  it  was  possible  to  determine  the  natural  fre¬ 
quencies  of  vibration  for  each  responding  mode.  These  natural  frequencies  then  can  be  compared  with 
those  predicted  by  the  NaTFREQ  computations  of  the  cable  with  attached  masses.  The  frequencies 
found  in  this  manner  actually  are  not  (he  undamped  natural  frequencies  but  the  damped  natural  fre¬ 
quencies.  However,  because  the  damping  is  so  low  in  air  the  two  are  essentially  the  same. 


From  free  vibration  decay  tests  logarithmic  decrements  were  estimated  from  the  acceleration-time 
histories  by  measuring  two  acceleration  amplitudes  separated  by  an  integer  number,  M,  of  complete 
cycles.  The  logarithmic  decrement  is  found  from  the  equation 


6 


A\ 


(19) 


where  A\  is  the  acceleration  amplitude  at  some  time  rt  and  is  the  amplitude  at  rw+ 1.  The  damp¬ 
ing  ratio  (  then  is  given  by 


49 


8 


(20) 


{ 


y/&2  +  (2  ir)2 


When  8  is  very  small  (8  «  1)  it  is  clear  that  £  =  8/2 n.  Details  of  the  analysis  procedure  are  given  by 
Moas  (11). 


The  natural  frequencies  of  several  modes  and  the  material  damping  were  measured  for  the  cable 
with  attached  masses.  The  damping  measurements  are  summarized  in  Table  4.  The  numbers  designat¬ 
ing  the  test  runs  in  the  table  correspond  to  the  numbers  in  Table  1  where  the  various  attached  mass 
configurations  are  given.  The  damping  often  was  found  to  be  nonlinear  over  the  range  of  in-air  dis¬ 
placement  amplitudes  (10).  Therefore,  when  two  damping  ratios  are  quoted  in  Table  4  for  a  particular 
test  run  the  first  value  corresponds  to  low  amplitudes  and  the  second  to  high  amplitudes  of  vibration. 
In  both  cases,  however,  the  material  damping  is  very  small  for  the  cable  with  up  to  seven  cylindrical 
attached  masses. 

The  measured  natural  frequencies  for  the  cable  with  attached  masses  are  listed  in  Table  5  for  the 
same  test  runs  listed  in  Table  4.  The  natural  frequencies  of  the  system  were  measured  by  exciting  the 
various  modes  in  both  the  vertical  (in-plane)  and  horizontal  (out-of-plane)  directions.  Generally  the 
results  are  comparable  in  both  cases.  A  discussion  of  the  data  reduction  sequence  employed  in  obtain¬ 
ing  the  natural  frequencies  is  given  in  Refs.  10  and  11.  It  was  found  in  general  that  the  vertical 
accelerometer  had  a  negligible  response  when  the  cable  was  excited  in  the  horizontal  plane,  and  vice 
versa.  The  natural  frequencies  were  determined  from  the  peaks  in  the  spectral  density  plots  of  the 
vibration.  Typically  there  were  no  mixed  mode  vibrations  and  only  one  peak  was  evident  in  the  spec¬ 
trum.  The  lowest  ( n  —  1)  mode  of  the  cable  apparently  was  influenced  by  the  sag  of  the  cable  (13). 
This  is  typical  of  the  vibrations  of  slack  cables,  where  the  influence  of  the  cable  sag  diminishes  as  the 
cable  vibrations  progress  to  higher  modes  (1,14).  Some  of  the  higher  mode  frequencies  were  probably 
influenced  to  some  extent  by  the  bending  stiffness  of  the  cable.  A  more  extensive  appraisal  of  cable 
sag  and  bending  stiffness  effects  should  be  undertaken  as  a  follow-up  to  the  present  study. 


A  typical  spectral  plot  is  given  in  Fig.  19.  This  measurement  was  by  the  accelerometer  pair 
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Table  4  —  Measurements  of  Cable  Material  Damping  (in  Air) 
1981  Castine  Bay  Field  Experiments, 

Marine  Cables  with  Attached  Masses: 

From  Refs.  10  and  13 


Test 

Number* 

Mode 

Direction  of 
Excitation 

Damping  Ratio,  £ 

1 

2 

Horizontal 

0.005 

3 

Horizontal 

0  001,  0.007 

3 

2 

Horizontal 

0  005 

3 

Horizontal 

0007 

3 

Vertical 

0  002.  0.003 

5 

2 

Horizontal 

0.004,  0.007 

2 

Vertical 

0  001.  0.003 

3 

Vertical 

0  001.  0.002 

7 

2 

Horizontal 

0  005.  0.006 

3 

Vertical 

t)  002.  0.005 

4 

Vertical 

0  002.  0.004 

9  . 

2 

Vertical 

0  004 

2 

Horizontal 

0  007 

3 

V'ertical 

0  002.  0  006 

3 

Horizontal 

0.004.  0.005 

11 

2 

Vertical 

0  008 

2 

Vertical 

0.005 

; 

2 

Horizontal 

0  008 

1 

3 

Horizontal 

0  008 

I 

3 

Vertical 

0  004 

1 

5 

Unknown 

0  004 

5 

Vertical 

0.003.  0.005 

19 

2 

Horizontal 

0  002.  0.003 

3 

Vertical 

0  001 

3 

Horizontal 

0  003.  0.004 

4 

Horizontal 

0  002 

M  roni  Tabic  I 
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Table  5  —  Measurements  of  Natural  Frequency  (In-Air) 
1981  Castine  Bay  Field  Experiments, 

Marine  Cables  with  Attached  Masses; 

From  Refs.  10  and  13 


Test 

Number* 

Mode 

Direction  of  NATFREQ-Predicted  Measured 

Excitation  Frequency,  Hz  Frequency,  Hz 

1 

2 

Horizontal 

1.72 

1.74 

i 

3 

Horizontal 

2.89 

3.036 

3 

2 

Horizontal 

1.59 

1.60 

3 

Horizontal 

2.36 

2.37 

3 

Vertical 

2.36 

2.40 

5 

2 

Horizontal 

1.73 

1.74 

2 

Vertical 

1.74 

3 

Vertical 

2.550 

2.63 

7 

2 

Horizontal 

1.51 

1.53 

3 

Vertical 

2.26 

2.34 

4 

Vertical 

2.98 

3.07 

9 

2 

Vertical 

1.57 

1.58 

2 

Horizontal 

1.61 

3 

Vertical 

3.00 

3  10 

i 

3 

Horizontal 

3.09 

i  , 

5 

Vertical 

4.04 

4.25 

u 

2 

Vertical 

1.60 

1.62 

2 

Horizontal 

1.62 

3 

Horizontal 

2.36 

2.41 

3 

Vertical 

2.46 

t 

; 

4 

Horizontal 

3.03 

3.18 

» 

5 

Vertical 

3.55 

3.80 

!  19 

2 

Vertical 

1.69 

1.75 

3 

Vertical 

2.62 

2.80 

t 

3 

Horizontal 

2.80 

i 

5 

Horizontal 

3.98 

4.30 

•From  Table  t. 
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located  at  the  position  L/ 4  along  the  cable.  The  light  ( L )  attached  masses  were  located  at  AT./6,  N  =■ 
1  to  5  For  this  example,  the  second  mode  (n  =  2)  with  five  light  masses  attached  to  the  cable,  the 
measured  natural  frequency  was  1  60  Hz.  The  natural  frequency  computed  with  NATFREQ  was  1.59 
Hz.  which  resulted  in  an  error  of  0.9  percent  (10).  The  error  between  the  measured  and  predicted 
natural  frequencies  ranged  from  this  value  to  about  11  percent.  However  in  most  of  the  cases  listed  in 
Table  5  the  error  between  the  measured  and  predicted  frequencies  was  less  than  5  percent 

The  measured  in-air  mode  shapes  for  the  cable  with  seven  light  attached  masses  (Test  Run  7)  are 
shown  in  Fig.  20.  The  normalized  measured  mode  shapes  are  denoted  by  the  distribution  of  individual 
displacement  spikes  at  the  seven  locations  of  the  accelerometer  pairs.  All  of  the  data  were  corrected  for 
rotation  of  the  accelerometers.  Superimposed  on  the  measured  displacements  is  an  equivalent  sine 
wave  The  measured  displacement  amplitudes  are  a  good  approximation  to  the  sine  wave  in  all  three 
modes  (n  —  2,  3  and  4).  Reference  to  Fig.  15(a)  readily  shows  that  the  NATFREQ-generatcd  mode 
shapes  also  closely  resemble  sine  waves  in  the  second  and  fourth  modes  where  the  masses  are  evenly 
distributed  over  a  half-wavelength  of  the  cable  vibration  pattern.  Even  for  the  third  mode,  where  the 
masses  are  not  evenly  distributed  over  the  half-wavelength,  there  appears  to  be  httle  vacation  from  the 
sinusoidal  distribution.  The  mode  shape  of  the  cable  does  not  resemble  a  nonsinusoidal  wave  form 
until  the  higher  modes  (n  *■  5  and  greater)  are  excited. 


E+07 


7.2  Strumming  of  Cables 


Laboratory-Scale  Experiments.  Experiments  were  conducted  by  Peltzer  (16)  to  study  the  effects  of 
shear  on  vortex  shedding  from  stationary  and  vibrating  marine  cables.  In  some  of  the  experiments,  a 
distribution  of  spheres  was  attached  to  the  cable  to  simulate  the  effect  of  attached  bodies  such  as  sensor 
housings  and  buoyancy  elements.  The  cable  of  11.4  mm  (.045  in)  diameter  had  a  length  to  diameter 
ratio  of  L/ D  -  107,  which  resulted  in  a  shear  flow  steepness  parameter  /§  =  0.0053  over  the  Reynolds 
number  range  Re  —  1.8  x  103  to  4  x  104.  Some  of  the  results  obtained  with  and  without  the  attached 
masses  are  summarized  here  for  comparison. 

The  spanwise  Strouhal  number  variation  .along  the  stationary  cable  in  the  shear  flow  is  shown  in 
Fig.  21.  There  are  ten  cells  of  constant  Strouhal  number  along  the  cable  span.  The  average  length  of 
the  cells  is  eleven  cable  diameters  and  the  change  in  the  Strouhal  number  from  cell  to  cell  is 

\StM  =  0.0086.  There  is  a  total  change  in  the  Strouhal  number  of  AS/W  -  0.067  across  the  span  from 

«#' 

-48  <  z/D  <48.  The  clarity  of  the  cellular  structure  is  due  to  a  careful  optimization  of  the  location  of 
the  hot-wire  probe  that  was  used  to  sense  the  frequency  of  vortex  shedding  from  the  cable  (15). 

The  cable  was  oscillated  in  its  first  mode  with  an  antinodal  displacement  amplitude  of  2  K  -  0.29D 
and  at  a  reduced  velocity  of  Vr  —  5.6.  The  vortex  shedding  pattern  in  Fig.  21  was  changed  by  the  oscil¬ 
lations  as  shown  in  Fig.  22.  The  vortex  shedding  was  locked-on  to  the  cable  vibration  over  the  central 
portion  of  the  cable  span  from  -14  <  z/D  <  30,  so  that  the  locked-on  cell  was  44  diameters  long. 
The  remainder  of  the  vortex  shedding  pattern  not  influenced  by  end  effects  was  stabilized  as  well.  Two 

A 

cells  were  increased  in  length  to  fourteen  diameters  each  and  no  fluctuations  in  the  cell  boundaries 
were  observed.  The  change  in  the  Strouhal  number  across  the  span  was  increased  to  A Stm  -  0.078. 

Ore  of  the  objectives  of  these  experiments  was  to  investigate  the  effects  of  attached  bluff  bodies 
on  the  vortex  shedding.  The  vortex  shedding  pattern  along  the  cable  with  five  spheres  (ping  pong 
balls)  attached  is  shown  in  Fig.  23.  The  flow  conditions  were  the  same  as  in  Figs.  21  and  22,  except 


57 


Fig  22  —  Strouhal  number  SlM  plotted  against  spanwise  distance  along  a  vibrating  flexible  cable  in  a  linear  shear  flow,  from 
Peluer  (16)  Reynolds  number  Re^  —  2.96  x  I03,  shear  flow  steepness  parameter  P  —  0  00S?  reduced  velocity  V,  -  5  6.  dis¬ 
placement  amplitude  (first  mode)  2  Y  —  0  290. 
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SPANWISE  DISTANCE,  z/D 


that  the  antinodal  displacement  amplitude  of  the  cable  was  2  Y  -  0.23Z). 

Several  points  concerning  the  vortex  shedding  from  the  cable-sphere  system  are  notable.  The 
three  central  spheres  (z/D  -  -20,  0,  20)  locked-on  to  a  submultiple  (one  half)  of  the  cable  vibration 
frequency.  The  lock-on  region  of  the  cable  strumming  in  Fig.  23  extended  from  -24  <  z/D  ^  29,  or 
over  53  diameters.  This  is  a  significant  increase  from  the  comparable  bare  cable  experiments.  For 
example,  when  2  Y  —  0.29 D  for  the  bare  cable  the  lock-on  region  extended  over  44  diameters,  and 
when  2  f*  0.23Dthe  lock-on  region  extended  over  34  diameters.  When  three  spheres  (z/D  =  -28, 
0,  +28)  were  attached  to  the  cable,  with  2  Y/D  =*  0.235,  the  lock-on  region  extended  over  61  cable 
diameters.  It  is  clear  that  the  addition  of  attached  bodies  along  the  cable  is  not  likely  to  deter  the 
resonant  cable  strumming  vibrations  even  in  a  shear  flow. 

Peltzer  (16)  observed  other  significant  features  relating  to  cable  strumming  in  a  shear  flow.  A 
maximum  separation  distance  of  twenty  cable  diameters  between  the  spheres  was  observed  that  would 
force  the  vortex  shedding  into  cells  of  constant  frequency  when  the  cable  was  stationary.  The  cellular 
structure  along  the  cable  with  spheres  attached  was  forced  into  a  pattern  that  was  appreciably  different 
from  the  bare  cable  in  the  shear  flow  when  the  spacing  between  the  spheres  was  twenty  diameters  or 
less. 

Towing  Channel  Experiments.  As  part  of  the  overall  NCEL  cable  dynamics  research  program  MAR 
Incorporated  was  funded  to  conduct  a  program  of  experiments  to  investigate  strumming  suppression 
and  the  effects  of  sensor  housings  (attached  masses)  on  the  overall  cable  response.  Two  reports  on  the 
results  obtained  recently  have  been  published  (2,3).  Some  of  the  results  obtained  are  summarized 
briefly  here  and  are  compared  with  previous  findings  from  the  NRL/DTNSRDC'/NCEL  cable  dynamics 
program.  The  tests  were  conducted  on  a  "strumming  rig"  at  the  DTNSRDC  that  was  employed  in  pre¬ 
vious  NCEL-sponsored  strumming  experiments  (1).  The  recent  MAR  Incorporated  experiments  and 
the  experimental  layout  are  described  in  detail  in  Ref.  3. 

Some  of  the  results  obtained  in  these  experiments  are  plotted  in  Fig.  24.  The  discrete  masses  in 
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all  cases  were  aluminum  sensor  housings  attached  at  various  locations  along  the  steel  cable  of  diameter 
D  **  0.175-in.  and  span  of  L  **  14.5-ft.  It  can  be  seen  from  the  results  in  Fig.  24  that  the  cable  was 
deployed  under  various  conditions  in  the  resonant,  cross  flow  strumming  regime  during  the  tests.  The 
results  plotted  in  the  figure  from  the  previous  DTNSRDC  experiments  were  obtained  with  a  bare  cable 
in  all  cases,  and  those  results  are  discussed  in  detail  in  Ref.  1 

The  attached  masses  did  not  deter  the  strumming,  but  instead  the  system  consisting  of  a  bare 
cable  plus  attached  masses  reached  higher  cross  flow  displacement  amplitudes  than  the  bare  cable  alone 
during  the  Mar  Incorporated  tests.  This  is  because  the  conditions  of  Mar’s  cable  experiment  were  at 
the  onset  of  the  resonant  strumming  regime  in  Fig.  24  while  the  attached  mass  experiments  extended 
well  into  the  resonant  region.  All  of  the  tests  were  conducted  in  the  range  of  cable  and  attached  mass 
properties  where  hydrodynamic  effects  are  dominant  (1)  and  even  the  addition  of  discrete  masses  has 
little  influence  on  the  large-displacement  cross  flow  strumming  effects.  All  of  the  frequency  spectra 
plotted  in  Ref.  3  give  clear  evidence  of  cross  flow  strumming  at  a  single  resonant  frequency. 

Cable  Strumming  in  the  Ocean.  Field  studies  of  the  strumming  behavior  of  marine  cables  were 
conducted  over  several  summers  of  the  Castine  Bay,  Maine  test  site  by  staff  members  and  students  of 
the  Ocean  Engineering  Department  at  MIT.  The  most  recent  experiments  are  described  in  this  report 
and  in  Refs.  8-12.  The  previous  experiments  are  described  in  Refs.  1.  5,  6  and  17. 

SEACON  II  was  a  major  undersea  construction  experiment  the  chief  goal  of  which  was  the  meas¬ 
urement  of  the  steady-state  response  of  a  complex  three-dimensional  cable  structure  to  ocean  currents. 
The  measured  array  responses  were  to  be  employed  in  a  validation  of  analytical  cable  design  models  and 
computer  codes  as  described  by  Kretschmer  et  al  (18).  A  second  goal  of  the  SEACON  II  experiment 
was  to  demonstrate  and  evaluate  new  developments  in  ocean  engineering  which  were  required  to 
design,  implant,  operate,  and  recover  fixed  undersea  cable  structures. 

The  SEACON  II  structure  consisted  of  a  delta-shaped  module  with  three  mooring  legs  It  was 
implanted  in  2900-ft  of  water  in  the  Santa  Monica  Basin  by  the  Naval  Civil  Engineering  Laboratory  dur- 
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Fig.  24  —  Strumming  displacement  amplitudes  for  a  small  diameter  cable  with  attached  masses:  data  from  Griffin  et  al  H)  and 
Kline  et  al.  (2).  Here  R  is  the  nsuhamof  the  in-line  and  cross-flow  strumming  components. 


STRUMMING  FREQUENCY,  f/Hz 


ing  1974  and  was  retrieved  during  1976.  The  top  of  the  cable  structure  was  positioned  450-ft  below  the 
water  surface  The  mooring  legs  were  4080-ft  long  and  each  arm  of  the  delta  was  1000-ft  long.  An 
artist’s  view  of  the  completed  structure  is  shown  in  Fig.  25.  Two  mooring  legs  were  attached  to  explo¬ 
sive  anchors  embedded  in  the  sea  floor  and  the  third  leg  was  attached  to  a  12500-lb  clump  anchor.  The 
entire  cable  structure  was  instrumented  in  order  to  collect  current  and  array  position  data. 

The  data  were  used  to  validate  the  NCEL  computer  code  DECELI.  This  code,  previously  called 
DESADE,  was  developed  at  NRL  by  Skop  and  Mark  in  1973.  The  delta  cables  experienced  uniform 
currents  over  their  respective  lengths  and  often  were  subject  to  cable  strumming.  These  strumming 
vibrations  led  to  increased  steady  drag  coefficients  and  static  deflections  as  described  in  the  next  section 
of  this  report  Details  of  the  SEACON  II  implantation,  design  and  recovery  are  given  by  Kretschmer  et 
al.  (18). 

Another  cable  strumming  experiment  (the  Bermude  Testspan)  was  conducted  by  the  U  S.  Navy 
from  December  1973  to  February  1974.  The  site  of  the  experiment  was  near  Argus  Island,  Bermuda. 
A  840-ft  long.  0.63-in.  diameter  electromechanical  cable  was  suspended  horizontally  in  the  water  at  a 
depth  of  92-ft.  The  cable  had  no  strumming  suppression  devices  attached,  but  it  had  numerous 
weights,  instrumentation  devices,  and  floats  distributed  over  its  length.  The  unlaired  cable  and  instru¬ 
mentation  were  similar  to  the  cables  which  made  up  the  horizontal  delta  segments  of  the  SEACON  II 
array.  Two  current  meters  were  suspended  near  the  mid-span  point  of  the  cable  Only  a  limited 
amount  of  data  were  obtained  from  this  experiment. 

Kennedy  and  Vandiver  (19)  have  analyzed  the  results  of  this  experiment  and  have  reached 
several  conclusions.  They  found  that  the  strumming  response  of  the  cable  was  typical  of  a  broadband 
random  process  and  that  resonant  and  nonresonant  lock-on  were  rare.  The  high  modal  density,  which 
resulted  in  responses  from  the  10th  to  the  150th  mode,  and  extreme  variations  in  current  speed  and 
direction  were  chiefly  responsible  for  the  broadband  response  of  the  test  span.  The  peak  rms  cross  flow 
displacement  amplitude  experienced  by  the  Bermude  Testspan  was  estimated  by  Kennedy  and  Vandiver 
to  be  Y  -  ±0.5D. 
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A  more  extensive  discussion  of  tnese  and  some  other  recent  cahle  strummins  field  experiments  is 


given  in  Refs.  1  and  17. 

7.3  Hydrodynamic  Drag 

An  important  side  effect  which  results  from  the  oscillations  of  structures  and  cables  due  to  vortex 
shedding  is  the  amplification  of  mean  drag  force  (or  the  drag  coefficient  CD).  The  drag  force 
amplification  measured  under  a  variety  of  conditions  has  been  summarized  in  a  recent  NCEL  report  (1) 
and  in  a  more  recent  and  related  paper  (17).  A  methodology  for  employing  these  measurements  in  the 
analysis  of  marine  cable  structures  was  developed  by  Skop,  Griffin  and  Ramberg  (20).  This  approach 
has  been  extended  to  the  case  of  flexible,  cylindrical  marine  structures  by  Griffin  (21,22).  A  step-by- 
step  method  for  approaching  this  problem  is  given  in  Refs.  1,21  and  22,  and  is  summarized  in  Appen¬ 
dix  C  of  this  report. 

The  measured  mean  drag  coefficients  ( CD )  for  several  strumming  cables  are  plotted  against  the 
Reynolds  number  (Re)  in  Fig.  26.  As  a  basis  for  comparison  the  typical  drag  coefficients  for  a  station¬ 
ary  circular  cylinder  and  several  nominally  stationary  braided  and  plaited  marine  cables  also  are  plotted 
in  the  figure.  The  relatively  large  scatter  in  the  stationary  cable  data  is  due  to  variations  in  the  low 
cable  tension  values  at  the  lowest  flow  speeds  (Reynolds  numbers).  All  of  the  cable  strumming  experi¬ 
ments  were  conducted  in  one  of  the  towing  channels  at  the  DTNSRDC.  The  dashed  line  was  simply 
faired  through  the  strumming  data.  It  is  clear  that  the  drag  coefficients  for  the  strumming  cables  are 
increased  substantially  (by  as  much  as  a  factor  of  two)  for  a  variety  of  Kevlar  cables  over  a  wide  range 
of  towing  speeds  or  Reynolds  numbers  between  Re  -  3003)  and  3G04). 

The  experimental  test  program  conducted  at  the  Castine  field  site  has  provided  a  more  extensive 
data  base  for  the  strumming  response  of  relatively  long  marine  cables  and  cylindrical  pipes  under  con¬ 
trolled  conditions.  The  drag  measurements  made  at  the  site  with  the  long  flexible  pipe  are  discussed  in 
a  separate  paper  by  Vandiver  (23).  A  time  history  of  the  cable  drag  coefficient  Cp  and  the  current  V 
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are  plotted  in  Fig.  27.  The  data  are  taken  from  a  typical  run  with  the  cable  only  and  not  attached 
masses.  From  the  approximately  3S  minutes  of  data  recorded  with  the  cable  strumming  it  is  clear  that 
the  hydrodynamic  drag  is  increased  consistently  from  the  level  that  is  typical  of  a  stationary  cylinder  or 
cable  ( CD  “  1  to  1.5).  The  measurements  of  CD  are  consistently  between  2  and  3  for  the  time  interval 
shown  when  the  current  velocity  is  near  1.2  kt.  The  strumming  response  of  the  cable  was  in  the  first 
six  (/i  -  l  to  6)  natural  modes. 

An  example  taken  from  one  of  the  more  complex  tests  (Run  20  of  Table  1)  is  shown  in  Fig.  28. 
Six  cylindrical  masses  were  attached  to  the  cable:  two  light  ones  at  LI 8  and  L/ 2;  and  four  heavy  ones 
at  L/ 3,  5Z./8,  3L/4  and  71/8.  The  RMS  strumming  response  data  shown  for  a  two  and  one  half  hour 
time  period  in  Fig.  28  were  recorded  at  31/4,  where  both  one  of  the  attached  masses  and  an  accelerom¬ 
eter  pair  were  located.  Several  segments  of  the  drag  and  strumming  displacement  amplitude  plots  are  at 
relatively  constant  levels.  These  probably  represent  time  periods  of  resonance  or  'lock-in'  between  the 
vortex  shedding  and  the  strumming.  Several  examples  taken  over  a  shorter  time  period  (448  seconds) 
are  now  discussed. 

A  segment  of  data  derived  from  Run  10  of  Table  1  is  shown  in  Fig.  29.  Two  heavy  masses  were 
attached  to  the  cable  at  Z./3  and  2L/3.  The  drag  coefficient,  the  current  speed,  and  the  horizontal  and 
vertical  rms  strumming  displacement  amplitudes  all  are  plotted  in  Fig.  29  for  a  typical  448  second  time 
period  during  this  test  run.  The  measurements  were  made  with  the  accelerometer  pa;r  at  5Z./8  along 
the  cable.  Frequency  spectra  for  each  of  the  horizontal  and  vertical  accelerometers  are  shown  in  Fig. 
30.  Three  stacked  graphs  are  shown  for  the  two  directions  and  each  graph  represents  a  strumming  fre¬ 
quency  spectrum  taken  from  a  two  and  one-half  minute  segment  of  the  data  file  for  this  test  run.  The 
displacement  output  derived  from  the  vertical  accelerometer  is  predominantly  at  a  frequency  near  3  Hz, 
with  a  smaller  spike  near  I  Hz.  The  spectra  derived  from  the  horizontal  accelerometer  are  character¬ 
ized  by  peaks  at  twice  those  of  the  vertical  accelerometer.  This  is  as  expected  since  the  fluctuating  drag 
force  in  the  horizontal  direction  has  a  frequency  of  twice  the  fluctuating  lift  foice  due  to  vortex  shed¬ 
ding.  The  predicted  natural  frequency  of  the  n  ■»  5  mode  for  this  test  run  is  3  1  Hz  It  is  likely  then 
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l-’ig  2?  —  A  lime  history  ol  the  drag  coefficient  Clt  and  ihe  current  velocity  V  recorded  during  the  1981  Casiine  Bay  field  test 
with  j  bare  cable  Cable  properties  as  rn  Table  2 


Fig.  28  —  A  two  and  one  half  hour  lime  record  of  the  hydrodynamic  drag  coefficient,  the  current,  and  the  vertical  and  horizontal 
strumming  displacement  amplitudes  for  a  cable  with  six  attached  masses  (Run  20  of  Table  1).  The  rms  strumming  amplitudes 
were  derived  from  the  accelerometer  pair  at  3L/4  along  the  cable. 


that  the  strumming  responses  shown  in  Figs.  29  and  30  are  either  at  or  near  this  resonant  frequency. 
The  response  spectra  also  may  contain  some  energy  at  the  Strouhal  frequency  of  the  flow  and  this  may 
account  for  the  twin  peaks  in  some  of  the  spectra 


The  coefficient  of  the  hydrodynamic  drag  force  was  amplified  considerably,  to  C/>  “  2.2  to  3.1,  as 
shown  in  Fig.  29.  This  is  a  substantial  increase  from  the  typical  drag  coefficient  of  Q>  —1  to  1.5  for  a 
stationary  cable.  The  vertical  displacement  amplitude  of  the  cable  strumming  varied  between  Y  — 
±0.2  and  0.6  in  (rms).  This  relatively  small  amplitude  at  the  5Z./8  location  is  consistent  with  vibra¬ 
tions  in  the  n  =  5  cable  mode,  since  a  node  of  the  vibration  pattern  then  is  to  be  expected  at  31/ 5 

nearby.  The  level  of  the  drag  amplification  would  appear  to  indicate  that  a  somewhat  larger  amplitude 

\ 

should  be  expected  at  the  antinodes  of  the  cabfe  vibration. 

Another  example  taken  from  one  of  the  more  complex  tests  (Run  14  of  Table  1)  is  shown  in  Fig. 
31.  Four  heavy  cylindrical  masses  were  attached  to  the  cable:  at  Z./8,  3Z./8,  51/8  and  7L/8.  The 
cable  tension  was  approximately  600  lb.  The  strumming  response  data  shown  for  the  448  sec  time 
period  in  Fig.  31  were  recorded  at  5Z./8  where  both  one  of  the  attached  masses  and  an  accelerometer 
pair  were  located.  The  vibration  level  over  the  time  of  the  test  run  was  approximately  Y  —  ±0.3  to 
0.45  in  (rms),  indicating  that  the  attached  mass  did  not  act  as  a  node  of  the  cable  system  vibration  pat¬ 
tern.  The  drag  coefficient  of  the  force  on  the  system  is  CD  —  2.4  to  3.2  which  rcprecenls  a  substantial 
amplification  from  the  stationary  bare  cable  value  of  Q>  =  1.2.  The  relative  contributions  of  the  cable 
and  the  masses  to  the  overall  drag  have  not  been  determined.  The  frequency  spectra  for  the  vertical 
and  horizontal  displacement-time  histories  in  Fig.  32  again  clearly  show  resonant  spectral  speaks  at  the 
vortex  shedding  frequency  (vertical)  and  twice  the  vortex  shedding  frequency  (horizontal)  The  cable 
strumming  vibrations  in  this  example  again  are  most  likely  to  be  in  the  fifth  mode  (n  —  5). 

\  third  example  is  given  in  Fig.  33.  For  this  case  (Run  16  of  Table  1)  six  heavy  masses  were 
attached  to  the  cable:  at  L/6,  L/3.  L/ 2,  5L/8,  3Z./4  and  1L/8  The  cable  strumming  response  meas¬ 
urements  were  made  at  L/6,  again  a  location  of  both  a  mass  and  a  pair  of  accelerometers.  The  drag 
coefficient  in  this  case  varied  from  about  CD  —  2.4  to  2  8  Over  most  of  the  time  history  shown  the 
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amplitudes  lor  the  cable  with  two  hca'-v  attached  cylindrical  masses  (Run  10  of  Table  I).  The  rms  displacement  amplitudes  were 
derived  from  the  outputs  of  the  accelerometer  pair  at  5L/&  along  the  cable 
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fig.  70(a)  —  Vertical  displacement  amplitude 

F  ig  30  —  Strumming  frequency  spectra  for  the  cable  with  two  attached  masses  (Run  10  of  table  0 
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|-ig.  31  —  A  lime  history  of  ihc  drag  coefficient.  :he  current  speed,  and  the  vertical  and  horizontal  strumming  displacement 
amplitude  for  ihc  cable  *uh  four  heavy  atiached  cylindrical  masses  (Run  14  of  Table  1)  The  rms  displacement  amplitudes  were 
derived  from  the  outputs  of  ihc  accelerometer  pair  ai  5£/8  along  the  cable 
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•ig  33  A  i true  hisr  ry  of  ihc  drag  inefficient.  the  curren*  speed,  and  the  vertical  and  horizontal  strumming  displacement 
mphtudes  for  the  cable  wuh  mx  heavy  attached  cylindrical  masses  (Run  16  ol  Table  l).  The  mis  displacement  amplitudes  were 
lerived  front  the  outputs  of  the  accelerometer  pair  at  /./ft  along  the  cable 


displacement  amplitude  varied  between  Y  -  0.1  and  0.4  in  (rms).  The  frequency  spectra  of  the  meas¬ 
ured  vertical  and  horizontal  strumming  displacement  amplitudes  shown  in  Fig.  34  are  non-resonant  with 
several  (five  to  six)  distinct  frequency  spikes.  It  is  likely  that  the  cable  tesponse  consists  of  several 
modes.  These  are  three  typical  examples  of  the  strumming  drag  for  the  cable  with  attached  masses. 
Similar  increases  in  drag  were  measured  for  alt  of  the  test  runs  with  masses  attached  to  the  cable. 

Several  test  runs  were  conducted  with  a  bare  cable  during  these  tests.  A  300  second  time  history 
for  one  such  test  is  shown  in  Fig.  35.  The  cable  was  resonantly  strumming  at  1.9  Hz  in  the  third  mode 
normal  to  the  current  and  nonresonantly  in  the  fifth  mode  in  line  with  the  flow  at  3.8  Hz.  The  vertical 
displacement  amplitude  is  approximately  Y  -  ±0.6  to  0.7  D(rms)  over  the  length  of  the  record  (7,8). 
The  average  drag  force  coefficient  on  the  cable  is  approximately  CD  =  3.2;  this  is  considerably  greater 
than  the  drag  coefficient  Q,  “  1  to  1.5  that  would  be  expected  if  the  cable  were  restrained  from  oscil¬ 
lating  under  these  same  flow  conditions.  The  drag  coefficint  on  the  strumming  cable  was  predicted  with 
the  equation  (see  Appendix  C) 

Co.  avg  “  Q,o  [1  +  1.043  (2  Trms/0)065). 

which  was  derived  from  the  original  results  of  Skop,  Griffin  and  Ramberg  (20)  The  strumming  drag 
coefficient  predicted  using  this  equation  is  in  the  range  CD  —  2.4  to  2.6  as  shown  in  Fig.  35.  This  is 
somewhat  below  the  drag  force  coefficient  measured  at  the  field  site,  but  the  predicted  values  are  rea¬ 
sonable.  The  predicted  drag  on  the  oscillating  cylindrical  pipe  was  virtually  indistinguishable  from  the 
measured  drag  coefficient  (8,23). 

The  N'ATFREQ  code  was  modified  as  described  in  Appendix  C  in  order  to  account  for  the  drag  on 
the  attached  masses.  This  now  is  done  in  the  two  subroutines  TTDRG  and  TWDRG  of  the  code.  The 
overall  mean  drag  coefficients  for  Runs  10,  14,  16  and  20  in  Table  1  were  computed  in  order  to  assess 
the  effect  of  the  attached  masses  on  the  overall  hydrodynamic  force.  The  final  results  are  given  in 
Table  6.  The  attached  masses,  between  two  and  six  in  number,  do  not  greatly  affect  the  overall  drag  on 
the  system  However,  the  stationary  cylinder  drag  coefficient  was  reduced  to  Cqo  “08  to  take  account 
of  the  small  aspect  ratio  (L/D  “  3.4)  of  the  cylindrical  masses.  For  a  cable  with  a  large  number  of 
attached  masses  there  is  likely  to  be  a  proportionately  greater  effect  on  the  drag.  The  comparison  is 
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Fig  M  -  Time  histories  of  the  measured  and  predicted  drig  coefficients  C0,  the  rms  cable  strummin«  displacement  amplitudes, 
and  the  current  speed  for  the  bare  cable  (7,8).  RMS  displacement  data  were  derived  from  the  accelerometer  pair  located  at  x  - 


only  qualitative  because  the  predominant  modes  of  response  corresponding  to  the  drag-time  histories  in 
Figs.  28,  29,  31  and  33  are  not  known.  The  average  drag  coefficients  for  the  mode  n  -  4  are  given  in 
the  table  as  examples.  These  are  representative  of  the  cable  strumming  responses  measured  during  the 
expei  iments.  The  computed  drag  coefficients  are  generally  within  or  close  to  the  range  measured  at  the 
Castine  Bay  test  site.  These  results  provide  a  reasonable  basis  for  concluding  that  the  NATFREQ 
computer  code  can  be  used  with  some  confidence  to  provide  engineering  estimates  of  the  mean  hydro- 
dynamic  drag  forces  on  marine  cables  with  arrays  of  discrete  masses  attached  to  them. 


Table  6 

Hydrodynamic  Drag  Forces 
1981  Castine  Bay  Field  Experiments, 
Marine  Cables  with  Attached  Masses 


Test 

Number 

Range* 

Drag  Coefficient, 
Cable  (Avg.)tt 

Drag  Coefficient, 

Cable  &  Masses  (Avg.)ttt 

10 

2.2-32 

3.00 

3.11 

14 

2.4-32 

2.76 

3  01 

16 

2. 1-2.9 

2.82 

3.13 

20 

2. 1-3.1 

2.83 

3.15 

+ Measured  at  (he  test  site. 

++Computed  using  NATFREQ;  Mode  number  n  —  Q>o  “  1-2- 
tf+Compuied  using  NATFREQ;  Mode  number  n  —  4-  Q>o  -  1-2. 
(cable),  CDo  m  0.8  (attached  masses) 


The  amplification  of  the  drag  measured  during  these  field  experiment;  is  comparable  to  the  drag 
increases  due  to  vortex  shedding  that  have  been  measured  in  recent  laboratory-scale  experiments  with 
cables  and  circular  cylinders  (23,24).  For  example,  Overvik  (24)  measured  drt.g  coefficients  of  CD  - 
2.5  when  a  model  riser  segment  (a  bare  circular  cylinder)  was  oscilla'ing  in  water  with  Y  -  ±1.1D. 
This  is  an  increase  of  230  percent  from  the  comparable  measurement  when  the  cylinder  was  restrained 
from  oscillating.  Griffin  et  al  (17)  have  summarized  recent  laboratory-scale  measurements  of  the  drag 
on  strumming  cables.  Comparable  increases  in  CD  were  commonly  obtained  in  the  tests  that  were 
reported. 
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7.4  ValkUtiM  of  NATFREQ 


The  results  obtained  in  air  can  be  used  with  the  in-water  cable  strumming  results  to  validate  the 
NATFREQ  computer  code  for  use  in  engineering  calculations.  T  natural  frequencies  for  the  cable  in 
water  were  derived  by  C.Y.  Liu  at  MIT  from  spectral  analysis  of  the  Castine  Bay  held  test  data.  In 
most  instances  the  measured  frequencies  plotted  in  figures  are  averages  derived  from  two  or  three  spec¬ 
tral  estimates. 

Measurements  in  Water.  Comparable  measurements  of  the  cable’s  natural  frequencies  to  those 
made  in  air  were  also  made  from  data  collected  while  the  cable  was  strumming  in  water.  The  natural 
frequencies  were  derived  from  power  spectral  density  plots  for  each  cable-added  mass  configuration  and 
were  computed  from  136  seconds  of  data.  Two  examples  of  the  spectral  density-frequency  plots  of  the 
vertical  acceleration  signal  are  shown  in  Figs.  36  and  37.  The  averaged  spectra  were  taken  over  a  136 
second  time  period,  or  over  4096  points.  The  vertical  lines  represent  NATFREQ- predicted  frequencies. 

The  accelerometer  pairs  located  at  1/8  and  L/2  were  employed  in  the  analysis.  It  was  anticipated 
that  the  LI 2  location  would  be  a  node  for  the  even  (/>  »  2,  4,  etc.)  cable  modes,  so  that  differences  in 
the  spectral  density  plots  from  the  two  locations  could  be  used  to  discriminate  between  odd  and  even 
cable  modes.  In  some  cases  this  proved  to  be  a  reasonable  approach.  Most  of  the  records  that  were 
analyzed  were  characterized  by  non-resonant  strumming  behavior,  so  that  each  frequency  spectrum  typ¬ 
ically  showed  a  distribution  of  peaks  from  which  several  participating  cable  modes  could  be  identified 
with  confidence. 

A  representative  comparison  of  the  computed  and  measured  natural  frequencies  is  given  in  Table 
7  for  those  test  runs  which  have  provided  reliable  accelerometer,  tension,  and  current  data.  For  the 
listed  cases  in  which  the  natural  frequencies  can  be  discriminated  reliably,  there  is  good  agreement 
between  the  predicted  and  measured  natural  frequencies.  Generally  the  differences  between  the  two 
sets  of  frequencies  are  less  than  ten  percent  in  all  cases,  and  the  horizontal  and  vertical  accelerometers 
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Kig  J6  —  Spectral  density  of  the  measured  natural  frequencies  in  water  for  the  Casline  Bay  cable  with  two  light  attached  masses 
(Run  2  of  Table  I).  The  vertical  lines  are  the  NATFREQ-predicted  frequencies. 
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fig  V  —  Spectral  density  of  the  measured  natural  frequencies  in  water  for  the  Castine  Bay  cable  with  seven  light  attached 
masses  (Run  8  of  Table  II  The  tertccal  lines  are  the  N  ATI  RfiQ- predicted  frequencies 


Table  7  —  Measurements  of  Natunfl  Frequency  (In-Water) 
1981  Castine  Bay  Field  Experiment, 

Marine  Cables  with  Attached  Masses 


Test 

Number* 

Mode 

N  ATFREQ-  Prednned 
Frequency,  Hr 

Measured 
Frequency,  Hz 

Vertical 

i  Horizontal 

2 

1 

0.63  " 

0.72 

2 

1.25 

1.40 

3 

2.10 

2.20 

4 

2.81 

2.70 

IRSI 

5 

3.12 

3.30 

3.27 

7 

4.53 

4.82 

4.76 

4 

2 

1.21 

1.33 

1.37 

3 

1.80 

1.85 

1.85 

4 

2.37 

2.60 

2.70 

7 

4.64 

4.83 

4.81 

8 

1 

0.56 

0.6 

0.6 

2 

1.13 

1.16 

1  16 

3 

1.67 

1.68 

1.69 

4 

2.22 

2.29 

2.37 

5 

2.73 

2.83 

2.87 

7 

3.56 

3.65 

3  67 

12 

1 

0.58 

0.76 

2 

1.16 

i  1.26 

1.25 

3 

1.71 

1.72 

1.72 

4 

2.22 

5 

2.62 

i  2.67 

1 

16 

2 

1.10 

1.10 

3 

1.62 

1.74 

1.74 

8 

5.30 

5.22 

•From  Table  l 
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at  the  two  measurement  locations  give  almost  identical  results  in  virtually  all  of  the  modal  determina¬ 
tions  listed  in  the  table. 

The  results  from  two  typical  test  runs  are  plotted  in  Figs.  38  and  39.  In  the  first  example  seven 
light  masses  were  attached  to  the  cable  at  mU 8,  m  -  1  to  7.  The  in-water  natural  frequencies  plotted 
in  Fig.  38  confirm  the  good  agreement  that  is  suggested  in  Table  7,  for  the  cable's  natural  modes  up  to 
n  —  7.  Another  example  is  given  in  Fig.  39  for  the  cable  with  two  light  masses  attached  at  L/ 3  and 
2JL/3.  The  measured  and  predicted  natural  frequencies  in  water  are  in  agreement  up  to  n  —  7  and  n  — 
II  (not  shown  in  the  figure).  The  latter  two  frequencies  were  identified  by  comparing  the  spectral  plots 
at  L/ 8  and  L/2 ,  and  by  noting  the  presence  of  peaks  in  the  spectra  from  the  vertical  and  horizontal 
accelerometers  at  L/2.  On  this  basis  it  reasonably  can  be  assumed  that  the  frequencies  are  odd  modes 
and  are  equal  to  n  ™  7  and  n  —  1 1  It  also  is  clear  from  Fig.  39  that  the  measured  natural  frequencies 
exhibit  the  same  discontinuous  behavior  that  is  predicted  when  the  number  of  cable  segments  is  equal 
to  the  mode  number  n. 

Several  in-air  natural  frequencies  also  are  plotted  in  Figs.  38  and  39.  As  had  been  noted 
earlier  there  also  is  good  agreement  between  the  measured  and  predicted  frequencies  in  these  cases. 
The  added  mass  effect  of  the  water  is  evident  from  the  computed  and  measured  results.  In  most  cases 
the  predicted  natural  frequencies  are  slightly  lower  than  the  measured  values.  An  added  mass 
coefficient  Cam  “  1  was  chosen  for  both  the  cable  and  the  cylindrical  masses.  This  value  may  be 
slightly  high  but  it  is  difficult  to  justify  another  choice  a  priori,  especially  since  the  measured  and 
predicted  frequencies  are  in  such  close  overall  agreement. 

It  is  clear  from  the  results  of  this  comparison  that  it  is  possible  to  predict  the  natural  frequencies 
of  a  taut  cable  with  attached  masses  to  better  than  ten  percent.  However,  at  the  higher  modes  the 
difference  in  frequency  between  modes  will  be  considerably  less  than  this  value.  Then  it  is  difficult  or 
impossible  to  positively  identify  which  spectral  frequency  is  associated  with  which  predicted  mode. 
Such  detailed  identification  would  require  more  extensive  knowledge  of  the  mode  shapes  of  the  system. 
As  a  practical  matter  this  requires  many  more  measurement  locations,  with  the  concomitant  expense  of 

additional  instrumentation  such  as  accelerometers  and  a  more  complex  cable  construction. 
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%  38  -  A  comparison  between  predicted  (NATFREQ)  and  measured  natural  frequencies  for  a  cable  with  seven  light  attached 
'.ndrical  masses  (Runs  ^  and  8  of  Table  II.  Measured  in  air.  0.  measured  in  water  •;  remaining  legend  as  in  Fig.  14 
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NATURAL  FREQUENCY,  fn/Hz 


10.0 


Fiit  M  —  A  comparison  between  predicted  (NATFREQ)  and  measured  natural  frequencies  for  a  cable  with  two  light  attached 
masses  (Runs  I  and  2  of  Table  I)  Measured  in  air.  0:  measured  in  water.#;  remaining  legend  as  in  Fig.  11. 
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8.  SUMMARY 


8.1  Findings  and  Conclusions 

Both  laboratory-scale  and  held  experiments  have  been  conducted  to  investigate  the  effects  of 
attached  masses  on  the  vortex-excited  strumming  response  of  taut  marine  cables.  A  comparison  has 
been  made  in  this  report  between  the  natural  frequencies  and  mode  shapes  computed  with  the 
NATFREQ  computer  code  and  measured  during  the  1981  Castine  Bay  field  experiment.  Also,  as  back¬ 
ground,  and  for  additional  comparisons,  summaries  are  given  here  of  the  NATFREQ  cable  analysis 
algorithm  and  other  recent  experiments  to  study  the  flow-induced  strumming  response  of  cables  with 
arrays  of  attached  masses. 

Twenty  test  runs  were  conducted  during  the  experimental  phase  of  the  program.  These  consisted 
of  ten  pairs  of  equivalent  tests  conducted  in  air  and  in  water.  The  measured  in-air  natural  frequencies 
are  in  good  agreement  with  the  NATFREQ  predictions  for  the  second  and  higher  (up  to  n  —  5)  cable 
modes.  The  first  mode  frequency  apparently  was  influenced  by  me  sag  of  the  cable.  The  measured 
mode  shapes  of  the  cable  vibrations  in  air  are  in  agreement  with  the  computed  mode  shapes,  but  only 
limited  mode  shape  comparisons  are  possible  due  to  the  existing  capabilities  of  the  code. 

Good  agreement  also  was  obtained  between  the  N  \T FREQ-predicted  natural  frequencies  and  the 
frequencies  measured  in  water.  It  is  clear  from  the  comparison  given  in  this  report  that  it  is  possible  to 
predict  the  natural  frequencies  of  a  taut  cable  with  attached  masses  to  better  than  ten  percent  How¬ 
ever,  at  high  cable  mode  numbers  the  difference  in  frequency  between  modes  often  is  considerably  less 
than  this  value.  Then  it  is  difficult  to  positively  identify  which  measured  cable  frequency  is  associated 
with  which  predicted  mode. 

The  results  available  to-date  tend  to  validate  the  NATFREQ  computer  code  as  a  reliable  engineer¬ 
ing  model  for  predicting  the  natural  frequencies  and  mode  shapes  of  taut  marine  cables  with  arrays  of 
discrete  masses  attached  to  them. 
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The  static  and  dynamic  analyses  of  a  marine  cabte  system  that  experiences  environmental  loading 
require  the  calculation  of  the  hydrodynamic  drag  forces.  Strumming  of  a  cable  segment  due  to  vortex 
shedding  increases  the  overall  mean  drag  force  and  results  in  a  corresponding  increase  in  the  cable  drag 
force  and  coefficient.  The  calculation  of  the  strumming  drag  requires  as  inputs  the  natural  frequencies 
and  mode  shapes  of  the  cable.  For  a  cable  segment  with  masses  attached  to  it  the  NATFREQ  code  pro¬ 
vides  this  information. 

Measurements  of  the  hydrodynamic  drag  forces  during  the  Castine  Bay  field  experiments  con¬ 
sistently  produced  drag  coefficients  in  the  range  CD  —  2  to  3.2  for  the  bare  cable  and  the  cable  with 
attached  masses.  This  is  a  substantial  amplification  of  the  drag  force  from  the  expected  level  for  a  sta¬ 
tionary  cable  (Co  —  1  to  1.5). 

The  drag  coefficient  on  the  strumming  cable  was  predicted  by  Vandiver  (23)  with  a  strumming 
drag  model  that  was  developed  at  NRL  as  part  of  the  overall  NCEL  cable  dynamics  research  program. 
The  predicted  drag  on  the  cable  was  within  about  20  percent  of  the  measured  drag.  For  a  complemen¬ 
tary  experiment  at  the  test  site  with  a  flexible  circular  steel  pipe,  the  predicted  drag  coefficients  were 
virtually  indistinguishable  from  the  measured  drag  coefficients  (23). 

The  NATFREQ  computer  code  contains  a  routine  for  computing  the  amplified  hydrodynamic  drag 
due  to  strumming-  This  routine  now  has  been  modified  to  account  for  the  drag  on  the  attached  masses 
(cylindrical  and  spherical).  Previously  only  the  strumming  drag  on  the  cable  segments  between  masses 
was  completed. 

8.2  Recommendation! 

Several  recommendations  have  been  developed  as  a  result  of  the  comparison  between  the 
NATFREQ  code  predictions  and  the  Castine  Bay  field  tests.  These  recommendations  are: 

*  Cable  sag  or  slack  cable  effects  often  play  an  important  role  in  the  dynamics  of  marine 
cables.  The  NATFREQ  code  is  limited  to  taut  cable  dynamics.  Consideration  should  be 
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given  to  the  development  of  a  comparable  code  for  the  calculation  of  the  natural  frequen¬ 
cies  and  mode  shapes  of  a  slack  cable  with  attached  masses.  At  least  one  such  code, 
called  SLAK,  presently  exists  in  rudimentary  form  (1). 


•  Many  of  the  data  records  from  the  1981  Castine  field  experiment  contain  lengthy  time 
segments  where  the  cable  strumming  is  nonresonant,  i.e.  the  vibrations  and  the  vortex 
shedding  are  not  locked-on  at  a  single  resonant  frequency.  Consideration  should  be  given 
to  developing  a  simple  but  still  effective  method  for  taking  nonresonant  strumming  effects 
into  account  in  determining  the  cable  response  and  the  strumming-induced  hydrodynamic 
drag  force  and  coefficient  Co- 
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Appendix  A 

A  LISTING  OF  THE  NATFREQ  COMPUTER  CODE 


oot-><">oor>oo<->r>or>oo<‘>c-»f>oOrir»<->*‘>o(‘)o<'>oooor»<^ 


FTN4,L 

PR  OCX  Ml  NEUNT 

C  PROGRAM  IEHNAT ( INPUT , OUTPUT  ,TAPE5=INPUT , TAPEA*OUTPUT ) 

C  TITLE!  •)  IS  81  COLUMNS  OF  TITLE  INFORMATION 

CDI  -KBUMK)  IF  N0NUNIF08H  CABLE  UEICHTS  AND  DIAMETERS  ARE  TO  BE 
SUPPLIES  FOR  EACH  SEGMENT 

1  .  IF  CABLE  HIEGHTS  AND  DIAMETERS  FOR  ALL  SEGMENTS  ARE 
SMC  AS  FIRST  SEGMENT 

CT  *0(ILAMK)  IF  NONUNIFORM  CABLE  TENSIONS  ARE  TO  BE  SUPPLIED  FOR 
EACH  CABLE  SEGMENT 

1  IF  CABLE  TENSION  FOR  ALL  SEGMENTS  IS  SAME  AS  FIRST 

SEGMENT 

NSEC  *  NO.  OF  SEGMENTS 
NPLOT  *  BLANK  -NO  PLOT  REQUIRED 
PLOT  -PLOT  OF  AMPLITUDE 
SCXL  -EACH  MODE,  SCRLL  FASHION 
PLEN  *  LENGTH  OF  DESIRED  AMPLITUDE  PLOT! IN) 

DEFAULT  *  21. 

CL!  1 )  s  LENGTH  OF  I-TH  CABLE  SEGMENT  (FT) 

D(I)  =  DIAMETER  OF  I-TH  CABLE  SEGMENT! IN) 

CDC(I)  =  DRAG  COEF  OF  I-TH  CABLE  SEGMENT 
UC(I)  -  TOTAL  HEIGHT  OF  I-TH  CABLE  SEGMENT  INCLUDING  ADDED  HATER 
HEIGHT  (LB/FT) 

TC(I)  *  TENSION  IN  I-TH  CABLE  SEGMENT  (LI) 

AH! I)  *  TOTAL  ADDED  HEIGHT  NT  END  OF  I-TH  CABLE  SEGMENT  INCLUDING 
ADDED  HATER  HEIGHT(LB).  AU(NSEG)  IS  ARBITRARY 
CDA(I)  a  CDSAREA  FOR  I-TH  ADDED  HEIGHT  IN  FTtt2 
STHT(I)  =  SPRUNG  HASS  AT  THE  END  OF  SPRING  AT  END  OF  THE  I-TH 
CABLE  SEGMENT  INCLUDING  ADDED  HATER  HEIGHT  (LB) 

MXMDS  =  MAX  NO  OF  NODES  TO  BE  FOUND  IN  SEARCH 
ACC  a  ACCURACY  PARAMETER  IN  MODE  SEARCH, USE  0.011  NOMINAL 
OMNAT  a  NATURAL  FREQUENCY  OF  THE  ATTACHED  SPRING-MASS  COMBINATION 
OMSTRT  a  BEGINNING  FREQ.  FOR  MODE  SEARCH  (RAD/SEC) 

OMSTOP  a  ENDING  FREQ.  FOR  MODE  S£ARCH(RAD/SEC) 

DELON  a  FREQ.  SEARCH  INCREMENT  (RAD/SEC) 

— CALCULATED  INTERNALLY  IF  DELOMaO 
UNITS  a  i  FREQUENCY  SEARCH  IN  RAD/SEC 
2  FREQUENCY  SEARCH  IN  HERTZ 
C  3  FREQUENCY  SEARCH  IN  FT/SEC 

C  4  FREQUENCY  SEARCH  IN  CM/SEC 

C  S  FREQUENCY  SEARCH  IN  KNOTS 

C  ISUPRS  a  |  FULL  PRINTOUT  OF  SEGMENT  RESPONSE 
C  l  SUPRESS  FULL  PRINTOUT 

C  ZETA  a  DAMPING  RATIO  FOR  CABLE  IN  AIR 

C  RHOU  a  HEIGHT  DENSITY  OF  FLUID  IN  LB/FTM3 

C  S(I)  «  ARC  LENGTH  TO  END  OF  I-TH  SEGMENT 

C  CP(I)  =  COUPLING  PARAMETER 

C  DU(I)  =  ARRAY  OF  DISTINCT  CABLE  DIAMETERS 

C  DS  a  DIA  OF  STRUMMING  CABLE  SEGMENTS 

C  STIF(I)aSPRING  STIFFNESS  AT  END  OF  THE  I-TH  CABLE  SEGMENT (cl /FT) 
COHHON/Ci/CL(  21). D<  21),HC(  21),TC<  2i),AH<  2i),NZX(  21>.«fSEC. 

1  ALP  (  21 ) 

C0MM0N/C2/A!  21), 8!  2l),PHAS(  21) ,PI ,BSIZE,AMPL(  21), RAMP  2D 
COMMON/C3/CD!  21) 
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C0HN0N/C4/BNAX,  JPLTNO 
C0HH0N/C6/STIF(  2i),STWT<  21) 

C0MH0N/C7/)N(  21),BNSIZE 

COHMQN/CU/CDC(  21),CDA(  21) >CP«  21)  ,DS,ZETA,ZETAE, RHOy,«mi, 

1CDAS<  21) 

DIKNSION  S(  21),DU<  21) 

INTEGER  TITLEUI), UNITS, CDI,CT 
REAL  KOMSR  T ,  KOHSTP , K DEL ON 
DOUBLE  PRECISION  0HCD,0HP,X1,X2,B1,B2,BP,IC 
C 

CStlttSET  FUNCTIONS  AND  CONSTANTS 
HERTZ(X)  »  X/(2.tPI) 

JPLTNO*! 

PI  =  ATAN2(I. ,-i. ) 

STRHN  =  1.21 
ND  =  1 

CtttttREAD/URITE  INPUT  DATA 

It  READ(S, ltll) (TITLE! J) ,J=1, 40) 

C  IF(EOF(S) >9999,11 

1111  FQRHAT(4IA2) 

C  11  CALL  DATE(IMY) 

C  CALL  TIHE(ITIN) 

UNITE (S,1II2)(TITLE(K),K-1, 41), IDAY,ITIN 

1112  FOR NAT ( 1H1 , 40A2/6H  DATE-, All, 7H  TIME-,AlO/iX,8l(iH=)//> 

ITIME  »  I 

IAPROX  a  « 

«EAD(S,1H3)CDI,CT 
1193  F0RNAT(2Ili) 

READ(S,1IB4)NSEG,NPLQT,PLEN 
1114  FORMAT(I10,A2,2X,F6.0) 

WRITE(6,108S)ND 

1915  F0RNAT(1X,12HDATA  SET  N0.,I3//17H  CABLE  PARAMETERS/) 

URITE(6,1886) 

1106  FORMAT (IX ,85HN0TE-  ALL  WEIGHTS  ARE  TOTAL  EFFECTIVE  WEIGHTS  AM)  INC 
1LUBE  EFFECTS  OF  ADDED  WATER  MASS) 

CtttltREAD  FIRST  CABLE  SEGMENT  INFO  HERE 

READ(S,1016)CL(l),D(l),CDC(l),UC(l),TC(t),AW(l),CDA(l).STIF(i), 

lSTWT(l) 

1016  FORHAT<F18  0,2F5.0,6F10.0> 

IF(CDI.EQ.i)URITE<6,l<07> 

1007  F0RHAT(7X,51HCABLE  WEIGHT  AND  DIAMETER  ARE  SAME  FOR  EACH  SEGFcNT) 
IFICT  EC  1)MRITE(6,10R8) 

1008  FORMAT (7X.38HCABLE  TENSION  IS  SANE  FOR  ALL  SEGMENTS) 

URITE(6.1009) 

1009  FORMAT!// , IX , 3HSEG ,5X ,6HLENGTH,3X, 19HARC  LENGTH, JX. 

1  4HD I AM , 3X , 4HDR AG , 3X , 9HWT/LENGTH , 3X , 7HTENS I ON , 

2  3X,8HATTACHED,3X,7HCDIAREA,3X,7HATT  SPR,4X,6HSPRUNC, 

3  3X,12H  SPRING-MASS/ 

4  2X,2HN0,6X,4H(FT) ,7X,4H(FT) ,6X,4H(IN) ,3X,4HCOEF,4X , 

5  7H(LB/FT)  ,5X,4H(LB)  ,5X,8HWT  (LB),3X,7H(FTM2).3X.7H(LB/FT) , 

6  3X , 8HNASS < LB ) , 3X , 1 1HN  FREQ  (HZ)/) 

COMMCHECK  HERE  FOR  NON  UNIFORM  WEIGHT, DIA  AND  TENSION 

IF(CDI  EQ  0  AND  CT  EO  I)  GO  TO  130 
CIOMICHECK  HERE  FOR  NON  UNIFORM  WEIGHT  AND  DIA  ONLY 
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IF(CBI.ER.O)  GO  TO  121 

cttstsaccx  hoe  for  nonuniform  tension  only 
IF(CT.EI.I)  GO  TO  ill 
CtttttUNIFORM  HEIGHT, II A  AM  TENSION 
M  ill  I  «  2, NSEC 

READ(S,10U)CL(I) ,AH(I),CDA(I) ,STIF(I) ,STHT(I) 

1111  FOSMT  (FIS .  •  ,30X ,4E1I .  I) 

UC(I)«HC(1) 

0(I)*D(1) 

CDC(I)  *  CK<1> 

TC(I)*TC(1> 
ill  CONTINUE 
GO  TO  141 

CMIWNONUNIFORH  TENSION  ONLY 
111  DO  111  I  =  2, NSEC 

READ(S,1912)CL(I),TC(I) ,AH(I),CDA(I) ,STIF(I) ,STHT(I) 

1112  F0RNAT(F19.I,2IX,5F1I.I) 

HC(I)=HC(1) 

CDC(I)<IC(1) 

D(I)*D<1) 

111  CONTINUE 
GO  TO  149 

CtttttNON  UNIFORM  HEIGHT  AND  DIA  ONLY 

129  DO  121  I  »  2, NSEC 

REAB(S,1I13)CL(I),D(I),CDC<I),MC<I) ,AU(I) ,CDA(I),STIF(I) ,STWT(I) 

1013  FORHATCFll  9,2F5.0,F10  0,10X,4F10 .0) 

TC(I)*TC(1) 

121  CONTINUE 
GO  TO  149 

CmttNON  UNIFORM  HEIGHT, DIA  AND  TENSION 

130  DO  131  I  *  2, NSEC 

READ(5,1I1&)CL(I) ,D(I),CDC(I),yC(I) ,TC< I) ,AH(I),CDA(I) , STIF (I ) , 
lSTUT(I) 

131  CONTINUE 

CtttttHRITE  OUT  CABLE  PARAMETERS 
141  CLT  ■  9  0 
AHT  *  II 
THIN  *  TC(l) 

UMAX  -  HC(1) 

SARC  *  II 
DHAX  =  10 
AH(NSEG)  <1.0 
DO  141  I  *  1,NSEC 
S(I)  =  01 
SARC  *  SARC+CKI) 

S<I>  =  SARC 
0MNAT  »  0 

IF  (STHT(I)  NE  D  OMNAT  =  S0RT<STIFa)t32.2/STHT(I))/<2.*PI) 
URITE(6,1014)I,CL(I),S(I>,D(I),CDC(I>,UC(I),TC(I),AU(I),CDA(I). 
1  STIF(I)  STHT(I)  OMNAT 

1014  F0RNAT(1X,I3,3X,F8  2,3X,F8.2,4X,F5.2)2X,F5.2)3X>F7.2,SX,F7.2, 

1  3X,F?  2,3X,F7  2,2X,F8  2,4X,F7  2,6X.F8  3) 

CttmCALCULATE  ALPHA(I) ,  CLT,  AHT.  THIN  .UMAX 
CLT  =  CLT  ♦  CL ( I > 
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MIT  •  MIT  ♦  MI<I)+STHT(I) 

IF  (TC(I)  LT.  THIN)  THIN  *  TC(I) 

IF  (UC(I)  CT.  UHAX)  UHAX  *  UC(I) 

ALP(I)=  SQRT(UC(I)/32.2/TC(I) ) 

IF  (Dd)  CT  DHAX)  DHAX  *  D(I) 

IF  (I  Efi  1)  WIN  -  WAX 
IF  (D(I)  LT.  WIN)  WIN  *  D(I) 

141  CONTINUE 

CtttttFIND  NO  OF  DISTINCT  DIA  AND  STONE 
X  *  1 

DU(1)»D(1) 

IFICDI.NE  0)  CO  TO  160 
DO  ISO  J  *  2, NSEC 
K  *  K  ♦  1 
DU(K>*  »( J) 

IL  *  J-l 
DO  14S  I  »  1, IL 

IF(ABS((D(J)-D(I))/D(J))GT  0015)  CO  TO  14S 
X  =  K  -  1 
CO  TO  ISO 
145  CONTINUE 
ISO  CONTINUE 
160  NOD  -  K 

IF  (MIT/CLT  CT.  UHAX)  UHAX  =  AUT/CLT 
NEAD<  S ,  1 0 15  )HXNDS  , OHSTRT ,  OHSTOP ,  DELOH ,  ACC .  WITS .  ISUPRS . 
1  ZETA,RHOU 

101S  FORHAT(I10,4F10  0,15,15, 2F10  0) 

IF(6HOM  EDO.)  RHOW  =  62  4 
Cttttt CHANCE  FROH  INPUT  SEARCH  UNITS  TO  RAD/SEC 

CALL  UNIT  (OHSTRT, OHSTOP, DELOH, UNITS, DHAX, DNIN) 
CSttttSET  DEFAULT  ACCURACY  TO  0  005  PERCENT 
IF  (ACC  LT  1.0E-7)  ACC  =  0.00105 
IF  (OHSTRT  EQ.  0.0)  OHSTRT  =  l  0E-10 
CtttttCALCULATE  DELOH  -  FREOUENCY  SEARCH  INCREHENT 
IF(D£LOH.NE. 0. >  CO  TO  210 
DELOH  =  SORT ( THINS32 . 2/UHAX ) /CLT /20  0 
210  URITE(6,1017)  NSEC 

1017  FORHAT  </24H  TOTAL  NO  OF  SECHENTS  =  ,13/) 

UNITE <6, 1030)  HXHDS 

1030  FORHAT  (1X,29HHAXIHUH  NO  OF  NODES  SOUGHT  =  ,13// 

1  15X,19H  NODE  SEARCH  LIHITS/UX,6H  LOWER, 4X,6H  UPPER/ 

2  11X,6H  L1HIT,4X,6H  LIHIT,2X,10H  INCREHENT) 

CtttttURITE  SEARCH  LIHITS  IN  RAD/SEC 

URITE(6,103i)  OHSTRT, OHSTOP, DELOH 

1031  FORHAT  (OH  RAD/SEC, 3F10  S> 

CtttttURITE  SEARCH  LIHITS  IN  HERTZ 

HOHSRT  =  HERTZ(OHSTRT) 

HOHFTP  *  HERTZ (OHSTOP) 

HDELOH  *  HERTZ (DELOH) 

URITE(6,1032>  HOHSRT .HOHSTP .HDELOH 

1032  FORHAT  (BH  HERTZ, 3F10  5) 

CtttttURITE  SEARCH  LIHITS  IN  FT/SEC 

FOHSNT  =  0HSTRTt(WIN/12.  )/(2.tPItSTRHN) 

FOHSTP  =  OHSTOP t(DHAX/l 2. )/(2.tPItSTRHN) 
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FBELOR  «  KL0H*(BHAX/12.)/(2.IPI«STIMN> 

URITE(4,1833)  FOHSRT  ,F0MSTP  ,FDELOH 

1133  FORMAT  (8N  FT/SEC, 3F10. 5) 

CtttttMRITE  SEARCH  LIMITS  IN  CM/SEC 

CONSRT  *  QNSTVTt(DMIN/12. )/(2.tPItSTRHN)t30.48 
COMSTP  =  «STOP*<DHAX/12.)/(2.OPISSTRHN>*30.48 
COELOM  >  BELOHt(DMAX/12. )/(2.tPI*STRHN)t30.48 
WRITE  (A, 1134)  COMSRT ,  COMSTP  ,CDELOH 

1134  FORMAT  <04  CM/SEC, 3F10. 5) 

CtmiWRITE  SEARCH  LIMITS  IN  KNOTS 

KOHSRT  *  QHSTRTXDHIN/12  )/(2.tPItSTRHN)>0 .5921 
KOMSTP  *  OMSTOPKDMAX/12. )/ (2. tPItSTRHN >11.5921 
K DELON  «  BEL0H0(DNAX/12.)/(2.tPItSTRHN)tQ  .5921 
MRITE(4,1035)K0H8RT,  KOMSTP  ,K  DELON 

1135  FORMAT  (W  KNOTS, 3F10.5) 

CtttttHRITE  INITIAL  CABUE  DAMPING  COEF.  AND  FLUID  DENSITY 
NRITE(6,1036)2ETA,RH0U 

1036  FORMAT!//, 32H  FRACTION  OF  CRITICAL  DAMPING  *  ,F7.4,//, 
1  17H  auID  DENSITY  »  ,lP£i0.4,10H  (LB/FTI3),/) 
CttttlBEGIN  SEARCH  FOR  NATURAL  FREQUENCIES 
OMP  =  -1.1 
NROOT  *  1 
ISERCH  x  R 
OHCD  =  OMSTRT 

311  IF  (OMCt  CE.  (OMSTOP  +  DELQM/iO. ))  GO  TO  SIR 
I FOUND  x  | 

311  CALL  SOLV(ONCD) 

BC  x  B(NSEC  ♦  1) 

IF(IFOU»  NE.I)  CO  TO  3S8 
IF (ISERCH. ME. I )  CO  TO  32R 
BP  «  BC 
ISERCH  x  i 

321  IF(BCtBP)340,340,330 
330  OMP  *  OHCD 
BP  x  |C 

OHCD  x  OHCD+  DELON 
GO  TO  301 
340  B1  x  BP 
B2  =  BC 
XI  =  OMP 
X2  x  (MCD 

3S0  IFOUND  x  IFOUND  +i 

IF  (IFOUND  GT  100)  GO  TO  400 
IF  (DABS(BC/BSIZE>  LE  ACC)  CO  TO  380 
IF(BltBC  LE  8  )  GO  TO  340 
B1  »  BC 
XI  *  OMCD 
GO  TO  370 
340  B2  *  BC 
X2  =  OMCD 

370  OMCD  *  Xi-U2-Xl)tBl/(B2-Bi) 

CO  TO  310 

CtttttNATURAL  FREQUENCY  FOUND 
3B0  CONTINUE 
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OUC  «  OHO 

CUtttCALCULATE  NOMINAL  AMPLITUDE  AND  PHASE  OF  EACH  SEGMENT 
CALL  AMPH 

CUtttCALCULATE  NO.  ZEROXINGS  AND  IDENTIFY  MODE  NO. 

CALL  ZEROX<OMCD,NNO> 

CtttttLOOP  FOR  ALL  DISTINCT  DIAMETERS 
DO  431  X  =  l,NOD 
DS  »  DUU) 

CtttttFORN  COUPLING  ARRAY 
DO  385  I  =  t ,NSEC 
CPCI)  =  0. 

IF(ABS((D(I)-DS)/DS).LT  8005)CP(I)  =  1. 

385  CONTINUE 

CtttttCALCULATE  IN  FACTOR  AND  OTHER  INTEGRALS 
CALL  INODE(OHC) 

CUtttCALCULATE  EFFECTIVE  DAMPING  FACTOR  AND  MASS 
CALL  DAMP 

CUtttCALCULATE  TRUE  RESPONSE  AMPLITUDE  FOR  EACH  SEGMENT 
CALL  RESMP 

CUtttCALCULATE  CD/CDO  FOR  EACH  CABLE  SEGMENT 
CALL  SCDRG(OHCD) 

CUtttCALCULATE  TOTAL  DRAG  COEFF  OF  CABLE  ALONE 
CALL  TTDRG  (CDT) 

CUtttCALCULATE  DRAG  COEFF  CORRECTION  FOR  STRUMMING  OF 
Cttttt  CYLINDRICAL  ATTACHED  MASSES 

CALL  TWDRC(DCDU) 

CUtttCALCULATE  TOTAL  DRAG  DUE  TO  STRUMMING  OF  CABLES 
Cttttt  AND  ATTACHED  MASSES 

CDUCT=CDT+DCDU 

CtttttHRITE  OUT  RESULTS  FOR  THIS  MODE 
HOHC  =  HERTZ(OMC) 

VELFT  =  0NCt<DS/12.)/<2.tPItSTRHN) 

VELCM  =  vaFT  t  31.48 
VELKTS  =  VELFT  t  0.5V21 
BN(MSEG+l>  *  BN(NSEG+1)/BNSIZE 
IF  (ITIME  .EQ.  0)  CO  TO  2111 
GO  TO  2112 
2011  HRITE(6,2I10) 

2011  FORMAT  (1H1) 

ITIME  =  1 
GO  TO  2013 

2012  WRITE(6,20Q1) 

2001  FORNAT(//lH  .8<15Htt - tt>/> 

2013  IF  < IAPROX  EQ.  0)  GO  TO  2014 
NRITE<6,2020> 

2020  FORMAT  (61H  NO  CONVERGENCE  THIS  MODE — APPROXIMATE  RESULTS  LISTEDt 

itttttt//) 

IAPRJX  *  0 

2014  URITE(&, 2015>MNO, OMC,HOMC, VELFT, K,VELCH,V€LKTS 

2015  FORMAT  <9H  MODE  NO  ,I4,/13H  FREQUENCY  *  ,F10.S.8M  RAD/SFC.3H  =  . 

1  F10.5,6H  HERTZ, 3X,27H  AVERAGE  STROUHAL  VELOCITY  ,F10  5, 

2  7H  FT/SEC/54X.21H  BASED  ON  DIAMETER  D(,I3.2H)  ,FiC  5.7H  CM/SEC/ 
380X,F10  5,6H  KNOTS) 

NRITE<6,2010)RMU,ZETAE 
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2114  FORMAT  (2W  EFFECTIVE  MSS  PARAMETER  =  ,1PE10.4, 

1  21H  EFFECTIVE  MMPINC  =  ,0PF10.6> 

IF(CDIEUi)  60  TO  3211 
60  TO  33R1 

3211  IF(lSUPRSECi)  GO  TO  421 
WRITER, 2002) 

2102  FQRNAT(/IX,7HSEGHEMT,3X,AHLENGTH,2X,10KARC  LENGTH, 2X, 

1  9HSEG  RESP,2X,8HUAV£  LNG,2X,UHN0  INTERNAL, 

2  2X.9HBRAC  COEF ,2X  ,8HATT  ACHED  ,3X ,9HMASS  RESP/ 

3  3X,2HN0,7X,4H(FT),AX,4H(FT),SX,9HANPL  (IN),5X,4H(FT> , 

4  6X , 6HZER0ES ,  12X ,4X ,5HNT  N0,5X,9HANPL  (IN)/) 

CO  TO  3S0I 

3311  MRITE(A,3305> 

3305  FORMAT  (108X,8H  SEGMENT/ 1X,7HSEGM£NT,3X,4HLENGTH,2X,10HARC  LENGTH, 

1  2X,  9HSEC  RESP,2X,8HUAVE  LNG,2X,liHN0  INTERNAL, 2X,9HDR AG  COEF, 

2  2X,8HATTACHEB,3X,9HMASS  RESP,7X,18H  STROUHAL  VELOCITY/ 

3  3X,2HN0,7X,4H(FT),AX,4H(FT),5X,9HAMPL  (IN), 

4  SX,4H(FT) ,&X,6HZER0ES,12X,4X,5HHT  N0,5X,9HANPL  (IN), 

5  3X,6HFT/SEC,4X,AHCN/SEC,5X,SHXN0TS/) 

3511  NSEG1  =  NSEC  -  1 

NZXT  =  I 
BMAX  *  B(i) 

DO  400  I  =  1,NSEG1 

IF (i(I)  GT  BMAX)  BMAX  -  B(l) 

IFCNPLOT  NE  2HPL  OR . NPLOT . EQ . 2HSC  )  GO  TO  35B0 
IF(NZX(1)  GT  I)  GO  TO  3571 
GO  TO  3580 

SEE  IF  CABLE  AMPLITUDE  EXCEEDS  NASS  AHPL  -  USE  GREATER 
VALUE  TO  SCALE  PLOT 

3571  SEC  *  II 

STEP  »  a(I)  /  (<NZX(I>  ♦1)17) 

3571  SEC  *  SEC  ♦  STEP 

IF(S£G  CT  CL(I))  GO  TO  3580 

AMP  «  RAMPd)  0  COS(SEC«ALP(I)tOMC-PHAS(I» 

IF  (AMP  GT  MAX)  BMAX  a  AMP 
CO  TO  3571 

3580  UVL  =  2  iPI/ALP(I)/OMC 
NZXT  =  NZXT  ♦  NZX(I) 

IF (CDI  EO  1 )  GO  TO  3791 
GO  TO  3801 

3701  HRITEU.,  2003)1,  CL(I),S(I>,  RAMPd),  NVL,NZX(I),CD(I),  I,  BN(I+1> 

2003  FORMAT(2X,I3,2X,F10.1,F10  1,2X,F10.4,FU.1,6X,I3,2X.F12  3, 

1  6X,I3,F134) 

CO  TO  400 

3801  VELFT  =  0MC»(D(I)/12  >/<2  0PI0STRHN) 

WELCH  a  VELFTI30.48 
VEUtTS  »  VELFTtO  .5921 

NRITE(A,  2100)1,  CL(I),S(I), RAMPd), MVL,NZX(I),CD(I).  I,  BN(I+1), 

1  VELFT,  WELCH,  VEUTS 

2100  FORMAT  (2X,I3,2X,Fl0.1,F10  i,2X,F10.4,Fii  1 ,4X,I3,2X,F12  3. 

1  4X,I3,F13  4,3<1X,F9  3)) 

400  CONTINUE 


100 


IFttlSEG  NE.i)  NZXT  >  NZXT  ♦  NZX(NSEC) 

WL  *  2.  *  P I/ALP ( NSEC )/OHC 
lF(CDI.EQ.l)  CO  TO  456 
GO  TO  478 

458  WRITE(6,2084)NSEG,CL(NSEG) ,S(NSEG) ,RAHP(NSEG) ,WWL,NZX(NSEG) 
1,C#<NSEC),BN(NSEC+1) 

2804  FORMAT  <2X,I3,2X,F18.1,F10.i,2X,Fi0.4,Fii.i,6X,I3,2X,F12.3,3X, 

*  8HBOUNDARY,Fll  4) 

CO  TO  498 

470  WELFT  =  0HCI(D(NSEG)/12. ) / (2 . *PI*STRHN) 

WELCH  =  WELFTI30 . 48 
WELKTS  *  WELFTIO .5921 

WRITE (6, 2151)  NSEC, CL (NSEC) ,S(NSEG) ,RAHP(NSEG) ,WWL,NZX(NSEG) , 

1  CD(NSEC) ,BN(NSEC+1 ) , WELFT, WELCH, WELKTS 
2158  FORHAT  <2X , I3,2X,Fi0  1 ,F10  1 ,2X ,F18 . 4 ,Fli . i ,6X , 13 , 2X ,F12  3 ,3X , 

1  8HB0UNDARY,F11.4,3(1X,F9  3)) 

490  UR ITE ( 6 , 2005 ) CLT , NZXT , CDT 

2885  FORHAT (/, IX, 5HT0TAL/  5X,F12  i,38X,I4,2X,F12  3,26H  (ADDED  NTS.  NOT 
1 INCLUDED)) 

495  NRITE(6,40Q0)  CDWCT 

4000  FORHAT ( 61X , F12 . 3 . 22H  (ADDED  UTS  INCLUDED)) 

GO  TO  430 

420  UR ITE (6,2030) CDT , BN ( NSEG+ 1 ) , IFOUND 
2030  FORHAT  (/2X,18H  AWERACE  CD/CD0  =  ,F10.3, 

1  26H  (ADDED  WTS  NOT  INCLUDED) ,5X,8H  AHPL  =  ,F12  S,5X, 

*  8H  ITER  =  ,110) 

430  CONTINUE 

CmttRESTART  SEARCH  FOR  NEW  NODES 
BP  =  B2 
OHP  »  X2 

OHCD  =  OHP  ♦  DELOH 

IFtNPLOT  E0.2HPL  )CALL  PLHPL  (TITLE, MNO,OHCD,WELKTS,PLEN,S> 
IF(NPLOT.E0.2HSC  )CALL  SCRLL  ( TITLE, HNO, OHCD, NROOT ,S) 

NROOT  =  NROOT  *  1 
IF ( NROOT. LE  HXHDS)  CO  TO  300 
C  SEARCH  COHPLETE  READ  NEW  DATA  SET 

S00  ND  =  ND  +  1 
ITIHE  =  0 
C  GO  TO  10 

GO  TO  9999 

CWMTERHINATE  SEARCH  IF  NO  CONWERGENCE 
600  IAPROX  *  1 
GO  TO  388 

C  CLOSE  PLOT  IFF  PLOT  WAS  PRODUCED 

9999  IF (  JPLTNO.NE.  0)  CALL  PL0T(i.,i  ,999) 

STOP 

END 

SUBROUTINE  UNIT  (Al,Bi,CfUNITS,Dl,D2) 

C0NH0N/C2/A(  21),B(  21),PHAS<  21>,PI,BSIZE,AHPL(  21),RAHP(  21) 
DIHENSION  0HEGA(3) ,DIAH<2) 

INTEGER  UNITS 

Cmt»SUBROUTINE  TO  CONWERT  FROH  INPUT  UNITS  TO  RAD/SEC  FOR  SEARCHING 
C  OHEGA(i)  =  OHSTRT 

C  0HEGA(2)  =  OHSTOP 
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C  0HECA(3)  *  Daw 
radix)  -  m  m 

STRHN  *  I  21 
OMECAd)  >  At 
0HECA(2)  *  B1 
0MECA(3>  «  C 
DIANd)  a  D1 
DIAH(2)  >  02 

GO  TO  (30. 40 ,51,61,71), WITS 
CRIMRUNITS  ARE  RAD/SEC 
31  RETURN 

CtttltUNITS  ARE  HERTZ 
4t  DO  42  Jsi ,3 
42  OHEGA(J)  a  RAD(ONECA(J)) 

CO  TO  til 

CmStUNITS  ARE  FT/SEC 
SR  DO  S2  1=1,2 

52  ONECA(I)  a  OM£GA(I)t2.tPItSTRHN/(DIAM(I>/12. ) 

OMEGA(3)  a  0MEGA(3)R2  tPItSTRHN/(DIAM(l)/12. ) 

GO  TO  111 

CtttltUNITS  ARE  CN/SEC 
60  DO  62  1=1,2 

62  OHEGA(I)  a  ONEGA ( I ) /3I . 48R2 . OPIOSTRHN/ ( DI AM( I )/ 12 . ) 

0HEGA(3)  a  ONEGA < 3 ) /30 . 4812 . <P I4STRHN/ (DIAH<1)/12.) 

GO  TO  111 

CMtftUNITS  ARE  KNOTS 
7fl  DO  72  I=i,2 

72  ONEGA(I)  =  QMEGA(I)/0  S92102 . OPIOSTRHN/ (DIAHt 1 )/12 . ) 

0HEGA<3)  a  ON£GA<3)/0 . 592142  0PI0STRHN/(DIAM( 1 )/12 . ) 

100  Ai  a  ONECA(i) 

Bi  a  0NEGA(2) 

C  =  ONEGA! 3) 

RETURN 

END 

SUBROUTINE  SOIV(ONC) 

COMMON/ Cl /CL (  2i),D(  21),HC(  21) ,TC(  21),AH(  21),NZX(  21). NSEC, 
1ALP(  21) 

COMNQN/C2/A(  21).B(  21>,PHAS(  21>,PI,BSIZE,AMPU  2i),RAMP(  21) 

COHHON/CS/AD,BD 

C0NH0N/C6/STIF<  21),STUT(  21) 

DOUBLE  PRECISION  AD(  2i),BD(  21),ARC,0HC,0MyT,C0EF 
AD(1)  =  1  000 
A(l)  a  A»(i) 

BD( 1 )s  l 
Bd)  a  Md) 

BSIZE  a  0  0 
DO  110  1  a  1.HSEG 
ARC  a  ALPd)IOMCtCL(I) 

WAVL  a  2  IPI/ALP(I)/OHC 
AN  *  DS0RT(AD(I)tt2+BD(I)tR2) 

IF  (<CL(1)  GT  NAVL/4  ). AND.  (AN. CT. BSIZE))  BSIZE  a  an 
IDU^l)  a  AD(I)*DSIN(ARG)+BD(I)IDCOS(ARC) 

1(141)  a  BDU41) 

IF(I  EO  NSEC)  CO  TO  100 
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ONHT  ■  It 

IF<STUTU )  GT  •  .010001  )0NMT*S8RT ( STIF<  1 )  132  2/STMT (I)) 
lF(DABS(ONMT/OMC>  GT. 1.9901  .0#  DABS(OHMT/OMC>.LT  09999) GO  TO  90 
HRITE<6,8I)0HC 

69  FORMAT <3X,4SHSPRUNG  MASS  ACTING  AS  ISOLATOR  AT  FREQUENCY  =  ,F10  5) 
90  COEF  =  STIF<I)/OMC/(i  B-(0MMT/0NC)902)-AW<I)00HC 

AD(I+1)  =  ALP<I)9<AD<I)9DC0S<ARG)-BD(I)9DSIN<ARC))9TC(I)/AIP<I+1)/ 
i  TC<I+1 )+CO£F  9BD<I+l)/ALP<I+i)/TC<I+i)/32.2 
A(Itl)  «  AD<I+i) 

IF  (tail)  LE  MAVL/4  )  AND  ( DABS( BD( I-fl > >  GT.  BSIZE)) 

1  BSIZE  =  DABS<BD(I+1)) 

100  CONTINUE 
RETURN 
END 

SUBROUTINE  AMPH 

COMMON/Cl/CK  21), D(  21),MC(  21>,TC(  2i),AM<  21),NZX<  21), NSEC, 
1ALP (  21) 

C0HM0N/C2/A<  21), B(  21),PHAS<  2i),PI,BSIZE,AMPL(  21),RAMP(  21) 

DO  191  I  *  1 ,NSEG 

Cl  =  SQRT(A( 1)992  ♦  B(I)992) 

C2  =  ATAN2(A(I),B(I)) 

AMPL(I)  =  Cl 
PHAS ( I )«  C2 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  IMODE(OMC) 

COMMON/Cl/CK  21) ,D(  21),WC(  21) ,TC<  21>,AM<  21),NZX<  21),N3EG, 
1ALP(  21) 

COMMON/C2/A<  21) ,B<  2l),PHAS(  21) ,P1 ,BSIZE,AMPL(  21),RAMP(  21) 
COMMON/C1 0/F2 ,F4,FG,F3,PIM 

COMMON/C 1 1 /CDC<  21>,CDA<  21>.CP<  21>,DS,ZETA,ZETAE,RHOW,RHU, 

1CDAS (  21) 

F4  =  0. 

F2  =  0 
F3  =  0 
FG  =  0 

DO  100  1=1  ,NSEG 

ARG  =  ALP ( I ) 90MC9CL ( I ) -PHAS ( I ) 

PAR  =  ALP(I)90MC9CL(I)/2.+<  SIN(2.9ARG>  ♦  SIN<2  9PHAS<I)>>/4 
F4  =  F4  ^ (AHPK I )9949(<C0S( ARG) 9939SIN<ARCHC0S<PHAS< I >)9939 
1  SIN(PHAS(I)))/4.+0  759PAR)/AIP(I)/0NC)9MC(I)+AN(I)9B<I+T)99 ' 

F2  =  F2  ♦  <AHPL<I)4029PAR/ALP<I)/0HC)9MC<I)+AM(I)4B(I*1)992 

FG  =  FG  ♦  CP ( I )9AMPL ( I ) 9929P AR /ALP ( I ) /ONC 

NPI  =  ALP  < I ) 90MC9CL  < I ) /P I 

EXTRAA  =  ALP<I)90MC9CL<  D-NPI9PI 

GAMMA  =  PHAS< I )*3.0Pl/2.-P I9AINT<<PHAS<  I )+3.9PI/2. )/PI) 

TERM1  =  SIN<EXTRAA-PHAS  < I ) )9<C0S(EXTRAA-PHAS< I ) )992+2  ) 

TERM2  =  SIN<PHAS< I ) J9<C0S<PHAS<  I ))992*2.) 

TERM3  =  <1  -CP < I ) ) ICDC < I ) 9D ( I ) 9 AMPL <1)993/36  / AIP(I  /OMC 

IF<EXTRAA  GT. GAMMA)  GO  TO  10 

F3  *  F3*TERM39<4  9NPI+ABS(TERM1+TERM2)) 

GO  TO  100 

10  F3  =  F3+TERM39  <  4  9NP I +ABS ( 2 . 9SIN<  GAMHA-PHAS  < I ) ) 9 
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i  (C0S(CA«MA  PHAS(I))Mg+2.)*TERH2-TERHl)) 

IIS  CONTINUE 
PIN  *  FVF2 
RETURN 
END 

SUBROUTINE  DNtf 

CONHON/Cl/CU  21), D<  21),«C<  21),TC(  21),A«<  2i),NZX(  21), NSEC, 
1ALP  <  21) 

C0NN0N/C2/A(  21), B(  21),PHAS(  21),PI,BSIZ£,ANPL<  2i),RANP<  21) 
CONMON/C1 1/F2 ,  F4 ,  FC ,  F3 ,  P  IN 

COWON/CU/CDC(  21),CDA<  2i),CP<  21) ,DS,ZETA,ZETAE,RHOW,RHU, 
1CDAS (  21) 

CtmtCALCULATE  EFFECTIVE  NASS  RATIO  PARANETER 
RNU  =  576 . SF2/P I /RH0U/DSII2/FG 
Ct  MISCALCULATE  ADDED  DAMPING  TERN 
ZADDED  *  I 
NSEC1  =  NSEC-l 
DO  21  I  -  1 ,NSEC1 

ZADDED  -  ZADDED+CDA<I)tABS(B(I+l)*«3> 

20  CONTINUE 

ZADDED  *  DStRHOUI(ZADDED+F3)/18 ,/PI/SQflT (F2)/SORT (F4) 
CMIMSEARCH  FOR  ZETAE 
ZETAE  =  ZETA 

25  VAR  =  (RNUIZETAE)tli  8 
PAR  -  1  *(9  6IVAR) 

DZETA  =  1  ♦  17.28  IZADDEDIRNUM1  8IZETAEII0  8/PARII2 

FZETA  =  ZETAE  -ZETA  -ZADDED/PAR 

ZETAEN=  ZETAE  -  FZETA/DZETA 

IFfABS(ZETAEN-ZETAE)  LT. . 1181)  GO  TO  30 

ZETAE  =  ZETAEN 

CO  TO  25 

30  ZETAE  =  ZETAEN 
RETURN 
END 

SUBROUTINE  RESNP 

COMNON/Cl/CU  21). D(  2i),MC(  21) ,TC(  21),AN(  2l),NZX<  2l ) , NSEG , 
1ALP  <  21) 

C0HN0N/C2/A(  21), B<  21),PHAS(  21) ,PI,BSIZE,ANPL<  2i),RANP(  21) 
CONNON/C?/BN<  21),BNSIZE 
C0NN0N/CII/F2 ,F4 ,FG ,F3 ,P IN 

COMMON/CU/CDC(  21),CDA(  21),CP<  21) ,DS, ZETA, ZETAE, RHOW, RNU. 
1CDAS<  21) 

ANAX  =  1  29/(1. 0*0  43«RHUtZETAE))l*3  35 
CONS  =  AHAXIDS/SQRTtPIN) 

DO  101  1=1, NSEC 
RAMP  < I )  =  CONStANPL(I) 

BN(I)  =  ABS(CONStB(I)) 

100  CONTINUE 

BNSIZE  *  ABS(CONSIBSIZE) 

BN ( NSEG* 1)  -  ABS(CONSIB(NSEG+i) ) 

RETURN 

END 

SUBROUTINE  SGDRG(ONC) 

CONMON/Cl/CK  21), D(  21),HC(  21),TC(  21),A«(  21),NZX(  21)  »SEC, 

104 


1ALP<  21) 

C0HH0N/C2/A(  21), B(  2i),PHAS<  21), PI ,BSIZE,AHPL(  21> ,RAMP<  21) 
CQHH0N/C3/CD(  21) 

COHMON/CU/CDC(  21)  ,CDA(  21)  ,CPt  21)  ,DS,ZETA,ZETAE,RHQH,RHU, 
1CDAS<  21) 

DOUBLE  PRECISION  QMC 
DO  111  I  =  1 ,NSEG 
R1  =  -1.  «  PHAS(I) 

PAR  =  ALP ( I )  I  ONC  *  CL (I ) 

R2  =  PAR  ♦  Rl  -  INT<PAR/PI)*PI 
CALL  SIHP(R1,R2,20,AREA) 

CD (I)  =  (1.  +  1 . 820241 (RAMP ( I ) /D < I ) ) *1 . 65* ( INT ( P AR/P I ) *2  2546+ 

1  AREA)/PAR)*CDC(I) 

100  CONTINUE 
RETURN 
END 

SUBROUTINE  SIHPtRl ,R2,N,AREA) 

F (X )=  ( ABS (  COS(X)) )I*0.6S 
DX  =  (R2-RD/2./N 
AREA  =  0 
X  =  R1 

DO  108  I  =  i,N 

AREA  =  AREA  »  (DX/3  )*(F(X)+4  *F(X+DX)+F(X+2.*DX> ) 

X  =  X  +  2  *DX 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  ZEROX(OHC,HNO) 

COHNON/C1/CL(  21), D<  2l),WC<  21).TC(  21)  ,AH(  2i),NZX(  21),NSEG, 
1ALP(  21) 

C0HH0N/C2/A(  21), B<  2i),PHAS(  21). PI  BSIZE,ANPL(  2»),RA«P(  ?1) 
C0NN0N/C5/AD,BD 

DOUBLE  PRECISION  ADI  2l),BD(  21) ,OHC.PR,PID,PHASD,PH 
PID  =  DABS (  DATN2I8  D0,-1  DO) ) 

NZX(i)  =  IDINT(ALP(i)*OHC*CL(i)/PID) 

MNO  =  NZX(l) 

NSEG1  =  NSEG-1 

IRNSEC.E0.2)  GO  TO  110 

DO  100  1=2 .NSEGl 

PH  =  ALP  < I ) tOMCICl ( I ) /P ID/2 . DO 

IF  (BD(I)*8D(I+i)  GT.  0.D0)  GO  TO  10 

NZX(I)  =  2»IDINT(PH>  +1 

GO  TO  SI 

10  NZX(I)  =  2*IDINT(PH+0  SD0) 

50  UNO  =  HNO  ♦  NZX(I) 

100  CONTINUE 

110  PHASD  =  DATN2 ( AD  <  NSEG ) , BD ( NSEG ) ) 

NL  =  ( 1+1DINTI DABS (PHASD) 12 . D8/PID) ) /2 
IF  (PHASD  .LT.  O  DD  NL  =  -t*NL 
PR  =  ALP ( NSEG )*OMC*CL (NSEG J-PHASD+l  SDO 
NR  =  ( 1+IDINT <DABS(PR )*2  DO/PID) )/2 
IF  (PR  LT  0  DO)  NR  =  -1*NR 
NZX(NSEG)  =  NL+NR 
HNO  =  NNO*NZX<WSEG> 
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RETURN 

END 

SUBROUTINE  TTNC  (CDT) 

CQNMON/Cl/CLt  21), D(  21),UC<  21>,TC(  21),AW(  2i),NZX(  21), NSEC, 
1RLP(  21) 

C0NM0N/C2/A(  21), 8<  21>,PHAS<  2i),PI,8SIZE,ANPl(  21),RAHP<  21) 
C0HM0N/C3/CD(  21) 

PR  *  0. 

CLT  =  I 

DO  111  I  =  1,NSEG 

PR  *  PR  ♦  Cl(I)*CD(I)t(D(I)/l2.l) 

CLT  =  CLT  ♦  CL(I)*(D<I)/12.8) 

111  CONTINUE 
CDT  *  PR/CLT 
RETURN 
END 

SUBROUTINE  TUDRC(DCDW) 

CONMON/Cl/CU  21), Ot  21 > ,WC<  21),TC<  21),AU<  21),NZX(  21),NSEG, 
iM.Pt  21) 

C0NN0N/C2/A(  21), B(  21),PHAS(  21) ,FI ; BSIZE , AHPLt  21).RANP(  21) 
CONMN/Cl/CDt  21) 

C0NN0N/C7/BNt  21),BNSIZE 

COHMON/Cit/CDCt  21) ,CDA(  21),CP(  21)  ,DS,ZETA,ZETAE,RHQW,RNU, 
lCDASt  21) 

NWSEC=NSEC-1 

PRD*I.I 

DO  81  1=1 ,NUSEC 

CDASt I ) =( t 1 . 0*1  82124ft BNt 1*1 )/D t I ) ) )tfO . 65)tCDA( I ) 
PRD=PRIKCDAStI) 

81  CONTINUE 
CLT*8.I 

DO  ill  1=1, NSEC 
CLT=CLT*CL(I)fD(I)/l2. 8 
108  CONTINUE 
DCDU=PRD/CLT 
RETURN 
END 

BLOCK  DATA 

COHNON/Cl/CLt  21), Dt  21),MC(  21) ,TC(  21),A«(  21),NZXt  21), NSEC, 
lALPt  21) 

CtWN0N/C2/A(  21), B<  2i),PHAS(  21) , PI, BSIZE, ANPLt  21),RAHP(  21) 
C0NN0N/C3/CD(  21) 

C0HH0N/C4/BNAX,  JPLTNO 
CONNON/CA/STIF(  21),STUTt  21) 

CONNON/C7/BNt  21),BNSIZE 

CONNON/Cll/CDC(  21),CDA<  21),CPt  21) ,DS,ZETA,ZETAE,RHOW,RNU, 
lCDASt  2D 

DOUBLE  PRECISION  ADt  2i),BD<  21) ,ARG,OHC,ONUT,COEF 

COMNON/CS/AD.BD 

COHNON  /C1I/F2.F4,FG,F3,PIH 

END 


Appendix  B 

INPUT  DATA  FOR  NATFREQ 
TEST  RUNS  1-20  FROM  TABLE  I 
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1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  1 


1  1 
3SCRL 

25. »0  1  25  12  7. 704E-01 

25. M 

25.10 

12  <M  50  0 


512  0  0  441E+01  3  500E-01 

0  441E+01  3  S80E-01 

0  441E+81  3  SI0E-01 

I  001  2  0 


.01 


0  0644 


1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  t  RUN  2 


1  1 

3SCRL 

25  80  1  25  1  2  0 

132E+01 

450  0  0  7S9E+01  3  SIOE-01 

25.00 

0  759E+01  3.500E-01 

25  10 

0  759E+01  3.5I0E-O1 

12  0  0 

SO  1 

0.001  2  0 

11 

64  0 

1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  3 


1 

6SCRL 


1 


12  50  1  25  1  2  ?  704E-01 

580.0  0  441E+01 

3  500E-81 

12  SI 

0  441E+01 

3  SOOE-Oi 

12.50 

0  441E+01 

3  5IOE-01 

12  51 

0  441E*01 

3  500E-01 

12  50 

0 . 441E+01 

3  SI0E-01 

12  50 

0  441E+81 

3.500E-01 

12  0  0  SO.O 

0  001 

2  0 

01 


0  0644 


1981  CASTINE  3AV  TESTS  *  NARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  4 
1  1 
fcSCRu 


12.50  1  25 

t  2  0  132E+01 

520  0  0  759E+01  3.500E-01 

12.50 

0  759E+01  3  500E-01 

12  50 

0  759E+01  3.5I0E-01 

12  SO 

0  759E+01  3  500E-01 

12.50 

0  759E+01  3.5I0E-01 

12  50 

0  7S9E+01  3  500E-01 

12  0 

0  SI  0 

0  001  2  0 

f  08 


1981  CAS7INE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  t  RUN  S 
i  1 

SSCRL 


9.375  1.2S  1  2  7.7I4E-01 

S32.0  0.441E+01  3.5I8E-01 

18.751 

0.441E+01  3.580E-Q1 

18  750 

0.441E+01  3.500E-01 

18  751 

0  441E+01  3.S80E-01 

9  375 

0  441E+01  3  500E-01 

12  0.0  SO.O 

0.001  2  0 

1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  t  RUN  6 
1  1 
SSCRL 


9.3751.25  1.2  0.132E+01 

475.0  0.759E+01  3.S00E-01 

18  750 

0  7S9E+01  3.S00E-01 

18.750 

0.7S9E+01  3  5O0E-O1 

18.750 

0.759E+01  3.S00E-01 

9.375 

0.759E+81  3.5I0E-01 

12  00  SO. 0 

0.001  2  0 

1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  7 
1  1 
8SCRL 


9.375  1.2S  1.2  7.704E-O1  584.0 

9  375 

9.375 

9.375 

9.375 

9  375 

9  375 

9.375 

12  0.0  50.0 


0  441E+01  3.S00E-01 
0  441E+01  3.5QQE-01 
0  441E+01  3.S00E-01 
0  441E+81  3  500E-01 
0.441E+01  3.5I0E-01 
0  441E+01  3  50OE-O1 
0.441E+01  3.5I0E-01 
0  441E+81  3.500E-01 

0.001  2  0  .01  0  0644 


1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSING'  *  RUN  8 
1  1 
SSCRL 


9.3751.25  1.2  0.132E+01  475.0 

9  375 

9.375 

9.375 

9.375 

9.375 
9  375 
9  375 

12  00  50  0 


0  7S9E+01  3.500E-01 
0.759E+S1  3.5Q0E-01 
0.759E+01  3.5I0E-01 
0  759E+01  3  500E-01 
0  759E+I1  3.5I0E-01 
0.759E+I1  3.S00E-  i 
0  759E+01  3  SOOE-fl 
O.7S9E+01  3  SQOE-^ 

0.001  2  01  M.O 
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1981  CASTINE  BAY  TESTS  >  MARINE  CABLES  WITH  ATTACHES  SENSOR  HOUSINGS  t  RUN  9 
1  1 
3SCRL 

25.10  1  25  1  2  7.7I4E-01  540.0  0.997E+01  3.S00E-01 

25.00  0.997E+01  3.580E-01 

25.00  8  997E+01  3  SO0E-O1 

12  0  0  SO  8  0.001  2  0  01  0  1644 


1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHES  SENSOR  HOUSINGS  *  RUN  10 

1  1 
3SCRL 

25.00125  1  2  0  132E+01 
25.00 

25  00 

12  80  50  0 


680  0  1.370E+81  3.5I0E-01 
1  370E+01  3  500E-01 
1  370E+01  3  50OE-01 

8  001  2  0  81  64  0 


1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  WITH  ATTACHES  SENSOR  HOUSINGS  *  RUN  11 


1  1 
6SCIL 


12.50  1  25  127  704E-01 

728.0  0  997E+01  3.5BOE-01 

12  SO 

0  997E+01  3  50OE-O1 

12  SO 

0  997E+01  3.500E-01 

12.  SO 

0.997E+01  3.500E-01 

12.50 

0.997E+01  3  5I0E-01 

12.50 

0  997E+01  3  SOOE-Bi 

12  0  8  58  0 

0  001  2  0 

1981  CASTINE  BAY  TESTS  •  MARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  12 


1  1 
6SCRL 


12.50  1  25 

1  2  0  132E+01 

580  0  1  370E+01  3  500E-01 

12.50 

1  370E*0i  3.500E-01 

12.50 

1  370E+01  3.500E-81 

12  SO 

1  370E+01  3  500E-81 

12.50 

1  370E*01  3.500E-01 

12  50 

1  370E+01  3  SOOE-Ol 

12  0 

0  50  0 

0  001  2  1 

1IO 


1981  CAST  IDE  BAY  TESTS  *  HARINE  CABLES  UITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  13 


i 

1 

5SCRL 

9.37S  1.25 

12  7.704E-01 

650.0  0  997E+01  3. SOOE-Oi 

18.750 

0  997E+01  3. SOOE-Oi 

18.750 

0.997E+01  3. SOOE-Oi 

18  7S0 

0.997E+O1  3.500E-01 

9  375 

0.997E+01  3  SOOE-Oi 

12  0 

0  50.0 

0.001  2  0 

1981  CASTINE  BAY  TESTS  1  MARINE  CABLES  UITH  ATTACHED  SENSOR 

i 

1 

5SCRL 

9.375  1.25 

1.2  0  132E+01 

650.0  i  370E+01  3. 500E-01 

18.750 

1.370E+01  3  500E-01 

18.750 

1  370E+01  3.5IOE-01 

18.750 

1.370E+01  3.500E-01 

9  375 

1  370E+01  3. SOOE-OI 

12  O 

.0  50  0 

Q.Ofli  2  0 

1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  UITH  ATTACHED  SENSOR 

1 

1 

7SCRL 

12.50  1  25 

1.2  7.704E-01 

732  0  0  997E+01  3.50QE-01 

12.50 

0  997E+01  3. SOOE-Oi 

12.  SO 

0.997E+01  3. SIOE-Oi 

9.375 

0.997E+I1  3. SOOE-Oi 

9.375 

0  997E+I1  3. SOOE-OI 

9.37S 

0.997E+01  3.500E-01 

9  375 

0.997E+01  3. SOOE-Oi 

12  0 

.0  50.0 

0.001  2  0 

0  0644 


t  RUN  14 


Si 


64.0 


01 


'  0644 


1981  CASTINE  BAY  TESTS  *  MARINE  CABLES  UITH  ATTACHED  SENSOR  HOUSINGS  »  RUN  16 


1  1 
7SCRL 

12  50  1.25  1.2  0 .132E+01  556. 

12.50 

12.50 

9.375 

9  375 

9.375 

9  375 


1.370E+01  3.500E->M 
1  378E+01  3  SOOE-OI 
1.370E+01  3.SIQE-01 
t.370E*0i  3  S00E-ft 
1  37OE+01  3  5I0E-I  S 
1.370E+01  3. SOOE-OI 
1  37IEMI  3  SI0E-G1 

0  001  2  0  01  64  0 


0.0 


1981  CASTINE  BAY  TESTS  t  NARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  17 


1  1 
6SCRL 

12.501.25  1.2  7.704E-01  800.0 

12.50 

12. 50 

12.50 

12.50 

12.50 

12  0  0  SO  .0 


0.997E+01  3. SOOE-OI 
0.441E+01  3. SOOE-Oi 
O.997E+01  3.500E-01 
0  441E+01  3 . SOOE-Oi 
0  997E+01  3.SI0E-O1 
0  997E+01  3  SOOE-Oi 

0.001  2  0  .01  0.0644 


1981  CASTINE  BAY  TESTS  *  NARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  18 


1  1 
6SCRL 

12.50  1.25  1.2  O  132E+01  765.0 

12.50 

12.50 

12.50 

12  50 

12.  SO 

12  00  50  0 


1.370E+01  3. 500E-01 
0  759E+01  3. SOOE-Oi 
1.370E+01  3  SOOE-Oi 
0  759E+01  3. SOOE-Oi 
1.370E+01  3.5I0E-01 
1.370E+01  3. SOOE-Oi 

0.001  2  0  .01  64.0 


1981  CASTINE  BAY  TESTS  *  NARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  19 


1  1 
7SCRL 

12.50125  1  2  7  704E-01 

12.  SO 

12.58 

9.375 

9  375 

9.375 

9  375 

12  0  0  SO  8 


800  0  0 .441E+01  3. SOOE-Oi 
0 .997E+01  3. SOOE-Oi 
0.441E+O1  3. SOOE-OI 
0.997E+01  3.58QE-01 
0.997E+01  3.5I0E-01 
0  997E+01  3.500E-01 
O.997E+01  3.500E-01 
0.001  2  0 


01 


0.0644 


1981  CASTINE  BAY  TESTS  *  NARINE  CABLES  WITH  ATTACHED  SENSOR  HOUSINGS  *  RUN  20 


1  1 
7SCRL 

12  SO  1  25  120  132E+01  800  fi 

12.50 

12.50 
9  375 
9  37S 

9.375 

9.375 

12  00  50  0 


0  7S9E+01  3  500E-01 
1  370E+01  3  500E-I1 
I  7S9E+01  3  5ME-01 
1  370E+01  3  500E-01 
1  370E+81  3.50OE-I1 
1  370E+I1  3  SOOE-Oi 
1  370E+01  3. 500E-01 

0  101  2  0  I.  64  0 


Appendix  C 

STRUMMING  DRAG  CALCULATIONS 


An  important  side  effect  which  accompanies  the  oscillations  of  structures  and  cables  due  to  vortex 
shedding  is  the  amplification  of  the  mean  hydrodynamic  drag  force  (or  the  corresponding  drag 
coefficient  C/>).  A  method  for  employing  these  measurements  in  the  analysis  of  marine  cable  struc¬ 
tures  was  developed  by  Skop,  Griffin  and  Ramberg  (201  '  This  procedure  has  been  improved  and 
extended  to  the  case  of  flexible,  cylindrical  marine  structures  by  Griffin  and  others  (1,21,22,25). 

A  program  of  tests  which  were  conducted  at  the  David  Taylor  Naval  Ship  R&D  Center  during  the 
1 940’s  and  for  which  the  results  were  released  a  few  years  ago  demonstrate  the  strong  resonance  due  to 
vortex  shedding  that  typically  takes  place  when  a  bare  circular  cylinder  or  cable  moves  in  steady  motion 
relative  to  the  surrounding  water  and  undergoes  large-amplitude  cross  flow  oscillations,  see  Ref.  21. 
This  cylinder  was  also  fitted  with  various  vortex  suppression  devices  in  order  to  investigate  their 
effectiveness  in  suppressing  the  cross  flow  oscillations.  The  drag  on  the  cylinuer  was  measured  over  a 
range  of  towing  speeds  up  to  10  knots  and  the  results  are  plotted  in  fig.  Cl.  A  clear  resonance 
occurred  near  1*4  knots  and  the  drag  force  (and  coefficient  Q>)  was  increased  by  a  factor  of  220  per¬ 
cent  at  a  towing  speed  of  4.25  knots.  At  this  and  nearby  towing  speeds,  the  cross  flow  displacement 
amplitude  of  the  cylinder  was  ±1.5  to  2  diameters.  When  the  cylinder  was  towed  at  speeds  above  and 
below  the  resonance,  the  usual 

Drag  a  (Flow  Speed)2 

dependence  was  obtained.  Methods  for  calculating  the  increased  mean  drag  forces  which  accompany 
the  vortex  shedding  are  described  in  the  discussion  that  follows.  Details  of  the  method  and  extensive 
comparisons  with  experimental  results  are  given  by  Griffin  *22).  and  by  Every,  King  and  Griffin  (25). 

The  drag  coefficient  CD  for  a  cable  or  structure  which  vibrates  due  to  vortex  shedding  is  increased 

'Numbers  in  parentheses  denote  references  listed  ai  the  end  of  this  reev  •: 
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60 


Fif.  Cl  -  Steady  drag  force  measured  on  a  free-ended  circular  cylinder  towed  steadily  through  still  water.  A  dear  resonance  in 
the  drag  due  to  the  vortex -excited  oscillations  normal  to  the  incident  flow  is  shown  between  relative  flow  speeds  of  V  -  3  and  6 
kts;  from  Griffin  121) 


as  shown  in  Figs.  26  and  Cl  of  this  report,  and  Fig.  2.9  of  Ref.  1.  The  ratio  of  CD  and  (the  latter 
is  the  drag  coefficient  for  a  cylinder,  cable  or  other  flexible,  bluff  structure  which  is  restrained  from 
oscillating)  is  a  function  of  the  displacement  amplitude  and  frequency  as  given  by  the  response  parame¬ 
ter  (20), 

w,  -  (1  +  2Y/D)(V,St)-'.  (Cl) 

Here  again  2Fis  the  double  amplitude  of  the  displacement,  V,  is  the  reduced  velocity  V/fD  and  5/ is 
the  Strouhal  number  /SZ>/  V.  The  ratio  of  the  drag  coefficients  is  given  by 

C0/ CO0  =1.  *S<1  (C2a) 

Co/ Coo  -  1  +  1.16U,  -  l)0  65,  w,  >  1  (C2b) 

which  is  a  least-squares  fit  to  the  data  in  Fig.  2.9  of  Ref.  1  The  equation 

?MAx/D  “  [1  +  0.43(27rS/2*s)]J35  (C3) 

from  Table  4.1  of  Ref.  1  can  be  combined  with  Eqs.  (C2a)  and  (C2b)  to  compute  the  unsteady  dis¬ 
placement  amplitudes,  the  drag  amplification  and  the  amplified  static  deflection  that  is  caused  by  the 
vortex  excited  oscillations.  The  local  displacement  amplitude  along  the  length  of  a  flexible  member 
such  as  a  cable  (in  the  »th  normal  mode)  is  given  by 

yiz)  *  Y,( r)  sin  <2irf,r) 

where 


F,(z)  -  Ff/7  wn>,(z). 


(C4) 


In  this  equation 


y, 


tr¬ 


ies) 


where  I ,  is  a  modal  scaling  factor  defined  in  Ref.  1.  These  equations  can  be  employed  to  iteratively 
compute  the  static  deflection  of  a  structure  or  cable  due  to  the  vortex-excited  drag  amplification  as 
shown  in  Refs.  1,  20  and  21. 


Every.  King  and  Griffin  (25)  recently  have  shown  by  comparison  between  sample  calculations  and 
the  experimental  data  reported  by  them  that  this  method  for  calculating  the  steady  drag  amplification 


and  deflections  is  quite  accurate  and  that  it  can  be  employed  with  some  confidence  to  evaluate  the 
hydrodynamic  loading,  deflections  and  material  stresses  on  marine  structures  and  cable  systems.  The 
and  unsteady  loads  and  deflections  must  be  considered  in  any  such  evaluation.  An  example  of  the 
comparison  reported  by  Every,  King  and  Griffin  is  given  in  Fig  C 2 

The  average  drag  coefficient  or  force  (with  respect  to  time)  can  be  derived  from  Eq.  (C2b).  The 
steps  required  to  do  this  were  carried  out  by  McGiothlin  (8)  who  failed  to  recognize  that  Eq.  (C2b)  is 

based  on  the  peak  amplitude  of  vibration  for  any  given  mode.  When  McGlothlin’s  derivation  is 

corrected  to  account  for  the  root-mean-square  vibration  amplitude  the  correct  result  is 

Co.avg  “  Cooil  +  1.043(2  WD)065! 

for  a  sinusoidal  mode  shape. 

The  hydrodynamic  force  due  to  cable  strumming  on  the  system  comprised  of  the  cable  and 
attached  masses  can  also  be  calculated  in  much  the  same  manner.  Then  the  total  drag  force  is  equal  to 
the  sum  of  the  drag  contributions  from  the  cable  strumming  and  from  the  attached  masses.  This  can 
be  stated  as 

Ft.  FCi  +  Fw,  (C6) 

Then  the  two  contributions  at  the  Ah  cable  segment  are 

Fc.,-  y  P  V'LACoc,  <C7«> 

and 

Fw.  ,  “  y  P  V1 LWIDWICDWJ  (C7b) 

for  each  of  a  total  of  N  segments  of  length  L,  and  diameter  D,.  There  are  N-  \  cylindrical  attached 
masses  of  length  LWl  and  diameter  DWi  situated  at  the  ends  of  the  individual  cable  segments.  Let  the 
total  drag  force  on  the  segment  be  defined  in  terms  of  the  local  cable  diameter  D,  and  segment  length 
Lt.  Then 

Ft.  ,  -  y  P  V'LACdt.  ,  (C8) 
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IN  LINE  DEFLECTION 


Fig  C2  —  The  steady  (ip  deflection  1,  of  a  cantilever  beam,  vibrating  due  to  vortex  shedding.  compared  with  the  measured 
values  for  a  cylinder  of  relative  density  SG  —  3  5  The  prediction  using  a  constant  CD  -  1.2  (stationary  cylinder!  i*  shown  by 
the  dashed  line,  from  Every.  King  and  Griffin  (25) 


»|7 


4. 


t 


and 


y  p  V*LiD,Cdt  ,  -  y  P  ^LD,Cdc  ,  +  y  p  K’LkAk  W 


(C9) 


When  this  is  simplified  and  the  summation  is  taken  over  Af  cable  segments, 
a-  ^  /v-i 

L,D,CDt  ,  «  LtDtCoc.  i  +  £wi  Oj^/  C0^ 

The  overall  average  drag  coefficient  Cor,  avg  can  be  defined  as 

N 


SA 


,  D, 


-or.  AVG- 


Then 


/  n  )_1  I N  N- 1 

Cot.  avg  “  2}  l,d\  hr  lac dc.  ,  +  J 


N  — ) 

Coc.  /  +  T  Lw,Dw,Cdw.  i 


(CIO) 


(Cl  2) 


(Cl  3) 


which  is  a  weighted-average  drag  coefficient  for  the  entire  system  based  upon  the  cross-section  areas  of 
the  cable  and  the  attached  masses  which  are  projected  into  the  incident  current.  For  a  spherical  mass 
the  cross-section  area  LAi  can  be  replaced  by  itD}/4  and  the  drag  coefficient  CDWi  by  the  appropriate 
drag  coefficient  for  a  sphere  when  an  appropriate  drag  amplification  algorithm  becomes  available.  Both 
drag  coefficients  Coc.  /  and  CDw.  ;  are  amplified  to  the  strumming.  And  the  cable  segment  vibration 
amplitude  and  the  vibration  amplitude  of  the  attached  mass  are  computed  in  the  NATFREQ  algorithm, 
so  that  the  drag  amplification  using  Eq.  (C2b)  can  be  computed  in  a  straightforward  manner.  It  is 
assumed  that  the  drag  on  the  attached  mass  (a  cylinder  of  L/D  —  3.4)  is  amplified  similarly  to  the  drag 
on  a  long  cable  or  cylinder.  The  summation  for  the  second  term  on  the  right  hand  side  is  taken  over 
N  -  1  segments  since  there  are  no  attached  masses  at  the  terminations. 


The  NATFREQ  code  was  modified  to  compute  the  drag  on  the  cable  system  in  this  way.  This 
now  is  done  in  the  two  subroutines  TTDRG  and  TWDRG.  The  drag  coefficients  computed  with  the 
modified  version  of  NATFREQ  are  compared  with  selected  examples  from  the  field  test  data  in  Sec.  7.3 
of  this  report. 
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